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Preface 


The  hydrodynamic  modeling  work  described  in  this  report  was  a  part  of 
the  New  York  Bight  Hydro-Environmental  Modeling  and  Monitoring  Study, 
authorized  under  Section  728  of  the  Water  Resources  Act  of  1986.  The  study 
was  performed  during  1989-1993  for  the  U.S.  Army  Engineer  District,  New 
York  (CENAN)  by  the  U.S.  Army  Engineer  Waterways  Experiment  Station 
(WES). 

Hydrodynamic  modeling  was  performed  by  personnel  of  the  Coastal 
Engineering  Research  Center  (CERC)  of  WES,  under  the  general  supervision 
of  Dr.  James  R.  Houston  and  Mr.  Charles  C.  Calhoun,  Jr.,  Director  and 
Assistant  Director,  respectively,  of  CERC.  Direct  supervision  was  provided 
by  Mr.  H.  Lee  Butler,  Chief  of  the  Research  Division  (RD),  CERC,  and 
overall  WES  Manager  of  the  New  York  Bight  project,  Mr.  Bruce  A. 

Ebersole,  Chief  of  the  Coastal  Processes  Branch,  RD,  CERC,  and  Dr.  Martin 
C.  Miller,  Chief  of  the  Oceanography  Branch,  RD,  CERC. 

This  report  was  prepared  during  1993  by  Dr.  S.  Rao  Vemulakonda  of 
CERC,  and  draws  heavily  on  an  earlier  user’s  guide  for  the  Chesapeake  Bay 
hydrodynamic  model  (Johnson  et  al.  1991b). 

CENAN  points  of  contact  during  the  study  were  Mr.  John  Tavalero, 

Ms.  Patricia  Barnwell-Pechko,  Ms.  Lynn  Bocamazo,  and  Mr.  Bryce 
Wisemiller.  Their  direction  and  assistance  were  invaluable.  Field  data  for  the 
1975-1976  period  were  collected  under  the  Marine  EcoSystem  Analysis 
(MESA)  program  and  were  made  available  under  contract  to  WES  by 
Dr.  Andrew  Stoddard  of  Creative  Enterprises  of  Northern  Virginia,  Inc.,  and 
Dr.  Gregory  Han  of  Han  &  Associates,  Inc. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was  Dr.  Robert 
W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


Conversion  Factors,  NON-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


Multiply 

By 

To  Obtain 

cubic  feet/sec 

0.02831685 

cubic  meters/sec 

feet 

0.3048 

meters 

1  Introduction 


The  New  York  Bight  (NY  Bight,  Figure  1)  is  generally  considered  to  be 
the  area  of  the  Atlantic  Ocean  extending  from  Cape  May,  New  Jersey  to 
Montauk  Point,  New  York.  Laterally,  the  NY  Bight  extends  about  160  km 
from  the  coastline  to  the  continental  shelf  break.  It  is  a  part  of  the  Middle 
Atlantic  Bight.  Depths  in  the  region  vary  from  3  m  near  Sandy  Hook,  New 
Jersey  to  900  m  in  the  Hudson  Canyon,  to  over  2,000  m  seaward  of  the 
continental  shelf.  Average  depths  are  approximately  60  m.  Circulation  of 
water  in  the  NY  Bight  is  influenced  by  astronomical  tides,  meteorological 
forcing,  flow  from  the  Hudson-Raritan  Estuaries,  density  gradients  due  to 
salinity  and  temperature  variations,  and  bathymetric  variations.  Large-scale 
oceanic  processes  affecting  the  entire  Middle  Atlantic  Bight  may  also  play  a 
role. 

The  NY  Bight  is  of  great  importance  for  a  variety  of  reasons,  including 
navigation,  recreation,  seafood,  and  disposal  of  wastes  and  dredged  materials. 
The  adjacent  coastal  states  are  heavily  populated  and  urbanized.  With 
increasing  population  and  industrialization  come  increasing  pressures  on  the 
ecosystem.  The  environmental  health  and  productivity  of  the  NY  Bight  and 
the  potential  impacts  of  natural  events  and  man-made  projects  on  the  Bight  are 
therefore  of  great  concern. 

This  report  is  a  product  of  the  NY  Bight  Study.  The  goal  of  this  feasibility 
study,  authorized  under  Section  728  of  the  Water  Resources  Act  of  1986,  was 
to  identify  a  means  of  developing  a  comprehensive  tool  for  effectively 
managing  the  resources  of  the  NY  Bight.  A  combined  hydrodynamic- 
environmental  modeling  technique,  used  together  with  a  monitoring  plan  and  a 
geographic  information  system,  was  identified  as  the  preferred  approach  to 
determine  potential  impacts  to  the  NY  Bight.  This  approach  was  adopted  and 
its  feasibility  was  demonstrated  during  the  course  of  the  study. 

Since  hydrodynamics  is  crucial  to  studying  other  aspects  such  as  water 
quality,  the  three-dimensional  hydrodynamic  model  CH3D-WES  (Curvilinear 
Hydrodynamics  in  Three  Dimensions)  was  selected  for  the  NY  Bight  Study. 
This  model  was  used  previously  by  the  Corps  of  Engineers  and  the 
Environmental  Protection  Agency  in  a  highly  successftil  joint  study  of  the 
Chesapeake  Bay  (Johnson  et  al.  1991a).  The  model  is  capable  of  using  a 
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Figure  1 .  The  New  York  Bight  study  area 


boundary-fitted  curvilinear  grid  in  the  horizontal  to  represent  complex 
geometries.  In  the  Chesapeake  Bay  application,  a  horizontal  (z-)  layer 
approach  was  used  in  the  vertical  direction.  Compared  to  Chesapeake  Bay, 
depths  in  the  NY  Bight  vary  considerably.  Therefore,  a  sigma-layer  version 
of  the  model  was  used.  In  this  approach,  the  same  number  of  layers  are  used 
in  the  vertical  direction,  in  both  deep  and  shallow  parts  of  the  study  area. 
Because  each  sigma  layer  represents  the  same  fraction  of  the  total  depth  at  a 
given  location,  sigma  layers  are  curved.  Even  though  the  main  interest  of  the 
study  was  in  the  Outer  Bight,  Long  Island  Sound  and  a  part  of  the  New  York 
Harbor  were  included  in  the  study  area  and  numerical  grid  for  the  sake  of 
completeness  and  to  represent  the  hydrodynamic  connections  between  the 
three  areas. 

The  numerical  grid  selected  for  the  study  (Figure  2)  has  high  resolution  in 
areas  of  greater  interest.  It  contains  a  maximum  of  76  cells  in  the  alongshore 
direction  and  45  cells  in  the  cross-shore  direction.  A  total  of  2,652  active 
horizontal  (water)  cells  are  used.  The  average  horizontal  grid  resolution  is 
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Figure  2.  New  York  Bight  computational  grid 

8  km  in  the  alongshore  direction  with  a  minimum  of  2  km  in  the  Hudson 
Canyon  and  a  maximum  of  17  km  near  Cape  May.  Average  resolution  in  the 
cross-shore  direction  is  6  km  with  a  minimum  of  200  m  in  the  East  River  and 
a  maximum  of  8  km  near  the  shelf  break.  The  Hudson  River  is  represented  . 
as  a  two-dimensional  laterally-averaged  body  of  water. 

The  hydrodynamic  model  was  calibrated  and  verified  for  April  and  May 
1976,  using  Marine  EcoSystem  Analysis  (MESA)  Program  field  data  for 
surface  elevations  and  currents.  A  6-month-long  extended  validation  was 
performed  for  the  period  April  1  -  October  1,  1976,  using  MESA  field  data 
for  salinity  and  temperature.  Results  from  the  simulation  were  used  to  drive 
successfully  the  Water  Quality  Model  for  the  NY  Bight.  A  demonstration  of 
the  hydrodynamic  model  was  performed  in  Long  Island  Sound  and  East  River 
for  a  72-day  period  starting  May  9,  1990,  using  the  National  Oceanic  and 
Atmospheric  Administration’s  (NOAA)  field  data  on  surface  elevations, 
currents,  salinity  and  temperature.  All  of  these  efforts  are  described  by 
Scheffiier  et  al.  (1994). 

The  present  report  is  a  user  guide  to  the  NY  Bight  hydrodynamic  model 
and  should  be  used  in  conjunction  with  the  report  by  Scheffner  et  al.  (1994), 
which  gives  more  details  of  the  theory  and  implementation  of  the  CH3D-WES 
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model,  as  well  as  numerical  simulations  performed  and  results  obtained.  The 
present  report  follows  the  following  outline.  After  this  introduction,  the 
CH3D-WES  model  is  briefly  described  in  Chapter  2.  Chapter  3  describes  the 
structure  of  the  model,  the  various  subroutines  used  and  the  function  of  each. 
In  Chapter  4,  the  setup  of  various  input  files  is  described,  using  an  application 
of  the  New  York  Bight  model  as  an  example.  Chapter  5  is  a  summary  of  the 

report.  Appendix  A  gives  model  parameters  and  dimensions  of  various  arrays 
used  and  Appendix  B  describes  the  format  for  the  primary  input  file. 
Appendixes  C-L  give  sample  input  files  for  the  application  given  in  Chapter  4. 
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2  CH3D-WES  Hydrodynamic 
Model 


The  numerical  hydrodynamic  model  CH3D-WES  (Curvilinear  Hydro¬ 
dynamics  in  Three  Dimensions  -  WES)  was  selected  to  provide  detailed 
hydrodynamic  flow  field  information  for  input  to  the  water  quality  or 
environmental  model.  For  convenience,  most  of  the  information  on  CH3D- 
WES  furnished  in  this  report  is  reproduced  from  Johnson  et  al.  (1991b).  The 
basic  model  (CH3D)  was  developed  by  Sheng  (1986)  for  WES,  but  was 
extensively  modified  later  in  its  application  to  the  Chesapeake  Bay  Study. 
These  modifications  have  consisted  of  implementing  different  basic  numerical 
formulations  of  the  governing  equations  as  well  as  substantial  recoding  of  the 
model  to  provide  more  efficient  computing.  As  its  name  implies,  CH3D-WES 
makes  hydrodynamic  computations  on  a  curvilinear  or  boundary-fitted 
planform  grid.  Physical  processes  impacting  circulation  and  vertical  mixing 
that  are  modeled  include  tides,  wind,  density  effects  (salinity  and  temperature), 
freshwater  inflows,  turbulence,  and  the  effect  of  the  earth’s  rotation. 

An  adequate  representation  of  the  vertical  turbulence  is  crucial  to  a 
successful  simulation  of  stratification  and  destratification.  A  second-order 
algebraic  turbulence  model  based  upon  the  assumption  of  local  equilibrium  of 
turbulence  is  employed.  The  boundary-fitted  coordinate  feature  of  the  model 
provides  grid  resolution  enhancement  necessary  to  adequately  represent  deep 
navigation  chaimels  and  irregular  shoreline  configurations  of  the  flow  system. 
The  curvilinear  grid  also  permits  adoption  of  accurate  and  economical  grid 
schematization  software.  The  solution  algorithm  employs  an  external  mode, 
consisting  of  vertically  averaged  equations,  which  provides  a  solution  for  the 
free  surface  displacement  and  unit  flow  rates  in  two  horizontal  coordinate 
directions  for  input  to  the  internal  mode,  which  solves  the  full  3D  equations. 


Governing  Equations 

The  governing  partial  differential  equations  are  based  on  the  following 
assumptions: 
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a.  The  hydrostatic  pressure  distribution  adequately  describes  the  vertical 
distribudon  of  fluid  pressure. 

b.  The  Boussinesq  approximation  is  appropriate. 

c.  The  eddy  viscosity  approach  adequately  describes  turbulent  mixing  in 

the  flow. 

The  basic  equations  for  an  incompressible  fluid  in  a  right-handed  Cartesian 
(x,y,z)  coordinate  system  (Johnson  et  al.  1991b)  are: 


du  dv  9w  „ 


(1) 


du  du^  duv 
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where 

(m,v,w) 

t 


f 

Q 

4) 

P 

p 

Ah,K„ 


g 

T 

S 


Po 


velocities  in  ix,y,z)  directions 
time 

Coriolis  parameter  defined  as  2Q  sin  (j) 
rotational  speed  of  the  earth 
latitude 

local  fluid  density 

local  fluid  pressure 

horizontal  turbulent  eddy  coefficients 

vertical  turbulent  eddy  coefficients 

gravitational  acceleration 

temperature 

salinity 

average  (reference)  fluid  density 


Equation  4  implies  that  vertical  accelerations  are  negligible  and  thus  the 
pressure  is  hydrostatic.  Equation  7,  commonly  known  as  the  equation  of  state, 
relates  the  local  density  to  the  local  temperature  and  salinity.  In  practice 
various  forms  of  the  equation  can  be  specified.  In  the  present  model,  the 
formulation  given  below  is  used: 


p  =  P/(a  +  0.698P) 


(8) 


where 

p  =  density  in  grams  per  cubic  centimeter 
P  =  5890  +  38r  -  0.375^  +  35 


a  =  1779.5  +  11.25r  -  0.07457'^  -  (3.8  +  0.017)5 


and  T  is  temperature  in  degrees  Celsius  and  5  is  salinity  in  parts  per  thousand 
(ppt). 

Non-Dimensionalization  of  Equations 

The  dimensionless  forms  of  the  governing  equations  are  used  to  facilitate 
comparisons  of  the  relative  magnitudes  of  various  terms  in  the  governing 
equations.  In  what  follows,  for  convenience,  the  same  symbols  as  before  are 
used  to  represent  dimensionless  variables. 


Chapter  2  CH3D-WES  Hydrodynamic  Model 


7 


External-Internal  Modes 


For  computational  efficiency,  the  solution  scheme  employs  a  "mode 
splitting"  technique.  In  this  procedure,  Equations  1-4  are  integrated  over 
depth  to  yield  a  set  of  vertically  integrated  equations  for  the  water  surface  ^ 

and  unit  flow  rates  U  and  V  in  the  x  and  y  directions  (external  mode).  Thus, 
the  time  step  for  the  rapidly  varying  external  mode  solution  is  limited  by  the 
surface  gravity  wave  speed  (the  Courant  condition).  The  major  purpose  of  the 
external  mode  is  to  solve  for  the  updated  water-surface  field  and  depth- 
integrated  flows  for  input  to  the  internal  mode  equations  from  which  the 
vertical  distributions  of  velocity,  salinity,  and  temperature  fields  are 
computed.  The  time  step  for  the  slowly  varying  internal  mode  solutions  is 
thereby  removed  from  the  Courant  condition  restriction  and  can  be  much 
larger  than  the  external  mode  time  step.  The  two  modes  together  provide  the 
full  3D  solution. 


Boundary-Fitted  Equations 

The  CH3D-WES  model  utilizes  a  boundary-fitted  or  generalized  curvilinear 
planform  grid  which  can  be  made  to  conform  to  flow  boundaries,  providing  a 
detailed  resolution  of  the  complex  horizontal  geometry  of  the  flow  system. 

This  necessitates  the  transformation  of  the  governing  equations  into  boundary- 
fitted  coordinates  (|,r;).  If  only  the  (x,y)  coordinates  are  transformed,  a 
system  of  equations  similar  to  those  solved  by  Johnson  (1980)  for  vertically 
averaged  flow  fields  is  obtained.  However,  in  the  CH3D-WES  model  not 
only  are  the  {x,y)  coordinates  transformed  into  the  (^,77)  curvilinear  system  but 
the  velocity  also  is  transformed  such  that  its  components  are  contravariant 
(i.e.,  perpendicular  to  the  (^,77)  coordinate  lines).  This  is  accomplished  by 
employing  the  definitions  below  for  the  components  of  the  Cartesian  velocity 
(n,v)  in  terms  of  contravariant  components  h  and  v  . 


U  =  XM  +  X  V 

V  =  y^u  +  yv 


(9) 


along  with  the  following  expressions  for  replacing  Cartesian  derivatives 


Chapter  2  CH3D-WES  Hydrodynamic  Model 


(10) 


/.  =  i  m',\  - 

/,  -  j  •  (ft^n 


where /is  an  arbitrary  variable  and  J  is  the  Jacobian  of  the  coordinate 
transformation  defined  as 


=  Vn  - 

With  the  governing  equations  written  in  terms  of  the  contravariant  components 
of  the  velocity,  boundary  conditions  can  be  prescribed  on  the  boundary-fitted 
grid  in  the  same  manner  as  on  a  Cartesian  grid  because  u  and  v  are 
perpendicular  to  the  curvilinear  cell  faces  (e.g.,  at  a  land  boundary,  either  u 
or  V  is  set  to  zero). 

The  vertical  dimension  is  represented  through  the  use  of  what  is  commonly 
referred  to  as  sigma-stretched  grid,  illustrated  in  Figure  3.  The  vertical  depth 
is  discretized  in  a  fixed  number  of  layers,  each  layer  equal  in  thickness  to  a 
fixed  percentage  of  the  local  depth. 

With  both  the  Cartesian  coordinates  and  the  Cartesian  velocity  transformed, 
boundary-fitted  equations  for  «  ,  v  ,  w,  5,  and  T  to  be  solved  in  each  vertical 
layer  are  obtained. 


Boundary  Conditions 


The  boundary  conditions  at  the  free  surface  are 


^du  dv^ 


dT 

IF 


TJ) 


dS 


=  0 


whereas  the  boundary  conditions  at  the  bottom  are 


(11) 
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Figure  3.  Sigma-stretched  grid 


dT 
dz  ’ 

where 


du 

3v 

Tz  ’ 

dz 

dS 

=  0 

dz 

'^STJ  ~ 

c  = 

W  = 


Pf.  = 
E.  = 

K  = 

T.  = 

'^bi^Tbr,  — 


=  (\  >  \)/p 


wind  stress  components 
surface  drag  coefficient 
wind  speed 

Vertical  Prandtl  Number 
Vertical  Ekman  Number 
surface  heat  exchange  coefficient 
equilibrium  temperature 
bottom  shear  stress  components 


(12) 


The  surface  drag  coefficient  is  computed  according  to  Garratt  (1977)  as 
follows 


C  =  (0.75  +  0.067  W)  X  10-3 


(13) 


with  the  maximum  allowable  value  being  0.003.  The  surface  heat  exchange 
coefficient,  K,  and  the  equilibrium  temperature,  T^,  are  computed  from 
meteorological  data  (wind  speed,  cloud  cover,  wet  and  dry  bulb  air 
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temperatures,  and  relative  humidity)  as  discussed  by  Edinger,  Brady,  and 
Geyer  (1974). 

Freshwater  inflow  and  water  temperature  are  prescribed  along  the  shoreline 
where  river  inflow  occurs,  however,  the  salinity  at  the  river  boundary  is 
specified  according  to  a  zero  spatial  gradient  assumption  (computed  from  the 
previous  time  step).  At  an  ocean  boundary,  the  water-surface  elevation  is 
prescribed  along  with  time-varying  vertical  distributions  of  salinity  and 
temperature.  Specified  values  of  salinity  and  temperature  are  employed  during 
flood  flow,  whereas  during  ebb,  interior  values  are  advected  out  of  the  grid. 
The  normal  component  of  the  velocity,  the  viscosity  and  diffusivity  are  set  to 
zero  along  solid  boundaries. 


Initial  Conditions 


When  initiating  a  run  of  CH3D-WES,  the  values  of  i^,  u,  v,  w,  U and  V are 
set  to  zero.  Values  of  salinity  and  temperature  are  read  from  input  files. 

These  initial  data  are  generated  from  prototype  measurements  at  a  limited 
number  of  locations.  Once  the  values  in  individual  cells  are  determined  by 
interpolating  from  the  field  data,  the  resulting  3D  field  is  smoothed. 

Generally,  the  salinity  and  temperature  fields  are  held  constant  for  the  first  few 
days  of  a  simulation  to  allow  the  flow  field  solution  to  spin  up. 


Numerical  Solution  Algorithm 


Finite  differences  are  used  to  replace  derivatives  in  the  governing  equations, 
resulting  in  a  system  of  linear  algebraic  equations  to  be  solved  in  both  the 
external  and  internal  modes.  The  external  mode  solution  consists  of  the 
surface  displacement  and  vertically  integrated  contravariant  unit  flows  UandV. 
All  terms  in  the  transformed  vertically  averaged  continuity  equation  are  treated 
implicitly  whereas  only  the  water-surface  slope  terms  in  the  transformed 
vertically  averaged  momentum  equations  are  treated  implicitly.  If  the  external 
mode  is  used  only  as  a  vertically  averaged  model,  the  bottom  friction  is  also 
treated  implicitly.  Those  terms  treated  implicitly  are  weighted  between  the 
new  and  old  time-steps.  Generally,  a  typical  value  of  0.55  for  the  time¬ 
weighting  parameter  yields  stable  and  accurate  solutions.  The  resulting  finite 
difference  equations  are  then  arranged  such  that  a  ^-sweep  followed  by  an  rj- 
sweep  of  the  horizontal  grid  yields  the  solution  at  the  new  time-step. 

The  internal  mode  consists  of  computations  for  the  three  velocity 
components  u,  v,  and  w,  salinity,  and  temperature.  The  only  terms  treated 
implicitly  are  the  vertical  diffusion  terms  in  all  equations  and  the  bottom 
friction  and  surface  slope  terms  in  the  momentum  equations.  Values  of  the 
water-surface  elevations  from  the  external  mode  are  used  to  evaluate  the 
surface  slope  terms  in  the  internal  mode  equations.  As  a  result,  the  extremely 
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restrictive  speed  of  a  free-surface  gravity  wave  is  removed  from  the  stability 
criteria  for  the  internal  mode  solution.  The  second  upwind  differencing 
scheme  of  Roache  (1976)  is  used  to  represent  the  convective  terms  in  the 
momentum  equations,  whereas  a  spatially  third-order  scheme  developed  by 
Leonard  (1979)  (called  QUICKEST)  is  used  to  represent  the  advective  terms 
in  the  transport  equations  for  salinity  and  temperature. 

It  should  be  noted  that  once  the  u  and  v  velocity  components  are 
computed,  they  are  slightly  adjusted  to  ensure  conservation  of  mass.  This  is 
accomplished  by  forcing  the  sum  of  u  over  the  vertical  to  be  the  vertically 
averaged  velocity  L7H  and  the  sum  of  v  over  the  vertical  to  equal  V/H,  where 
H  is  the  total  water  depth. 


Turbulence  Parameterization 


Vertical  turbulence  is  handled  by  using  the  concept  of  eddy  viscosity  and 
diffusivity  to  represent  the  velocity  and  density  correlation  terms  that  arise 
from  a  time  averaging  of  the  governing  equations.  These  eddy  coefficients 
are  computed  from  mean  flow  characteristics  using  a  simplified  second-order 
closure  model  originally  developed  by  Donaldson  (1973).  The  closure  model 
has  been  further  developed  and  applied  to  various  types  of  flows  by  Lewellen 
(1977)  and  Sheng  (1982,  1986,  1990).  The  procedure  assumes  local 
equilibrium  of  turbulence.  For  more  details,  the  interested  reader  should  refer 
to  these  references  and  to  Johnson  et  al.  (1991a). 


Computational  Grid 

A  staggered  grid  (Figure  4)  is  used  in  both  the  horizontal  and  vertical 
directions  of  the  computational  domain.  In  the  horizontal  direction,  a  unit  cell 
consists  of  a  ^point  in  the  center  (^.^.),  a  U-point  to  its  "west"  iU,j),  and  a  V- 
point  to  its  "south"  (V-j).  In  the  vertical  direction,  the  vertical  velocities  are 
computed  at  the  "full"  grid  points.  Horizontal  velocities,  temperature, 
salinity,  and  density  are  computed  at  the  "half  grid  points  (half  vertical  grid 
spacing  below  the  full  points). 

Two  arrays,  NS  and  MS  are  set  automatically  by  the  model  and  used  to 
flag  the  grid  cells  in  the  horizontal.  The  array  NS  indicates  the  condition  of 
the  "west"  and  "east"  cell  boundaries,  whereas  the  array  MS  denotes  the 
condition  of  the  "north"  and  "south"  cell  boundaries  (Figure  5). 

The  boundary-fitted  grid  shown  in  Figure  2  was  developed  to  provide  a 
high  resolution  representation  of  the  complex  geometry  of  the  coupled  New 
York  Bight,  New  York  Harbor,  and  Long  Island  Sound  model  domain.  The 
computational  grid  contains  a  maximum  of  76  cells  in  the  alongshore  direction 
and  45  cells  in  the  offshore  direction.  There  are  2,652  active  horizontal  grid 
cells  and  10  vertical  layers,  resulting  in  26,520  computational  cells.  The 
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Figure  4.  Staggered  grid 

Hudson  river  is  parameterized  as  a  two-dimensional,  laterally-integrated  body 
of  water  and  modeled  as  a  river  boundary  input  to  the  Bight. 


Depths,  specified  at  the  cell  comers,  were  extracted  from  bathymetric  data 
gathered  from  the  National  Oceanographic  and  Atmospheric  Administration’s 
(NOAA)  nautical  charts  based  on  Mean  Low  Water.  Although  the  resolution  in 
Long  Island  Sound,  New  York-New  Jersey  Harbor  and  Hudson-Raritan  estuary 
system  is  coarse  for  detailed  study,  the  overall  grid  resolution  is  sufficient  to 
provide  insight  into  the  circulation  and  transport  processes  in  the  Bight. 
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3  Structure  of  the  New  York 
Bight  3D  Hydrodynamic 
Model 


The  New  York  Bight  model  has  a  main  program  as  well  as  several 
subroutines.  Subroutines  governing  model  setup  are  called  from  the  main 
program  while  subroutines  governing  computations  are  called  from  subroutine 
CH3DM2.  Each  of  these  is  listed  below  with  a  description  of  its  function. 
Entry  points  in  subroutines  are  also  noted.  Note  that  the  UNIX  operating 
system  used  on  the  Cray  Y-MP  computer  and  the  Silicon  Graphics’  IRIS 
Indigo  Workstation,  where  the  CH3D-WES  model  is  run,  is  sensitive  to  lower 
and  upper  case  letters  in  file  names,  commands,  etc.  Two  INCLUDE  files, 
nyb.inc  and  ch3d.inc,  are  needed.  They  are  used  to  set  up  parameters, 
dimensions  of  various  arrays,  and  COMMON  blocks.  During  model 
compilation,  these  files  are  inserted  wherever  the  INCLUDE  statements  in  the 
source  code  call  for  the  files.  These  files  are  discussed  in  detail  in  Chapter  4, 
along  with  input  files. 


CH3D  The  main  program. 

CH3DIR  Reads  data  from  main  input  file,  FILE  4  (see  Appendix  B), 
which  controls  computations,  input,  and  output.  Various 
constants  are  computed,  and  the  vertical  (a-)  layer 
thicknesses  are  set. 

CH3DTR  Reads  (jc,y)  coordinates  (ft)  and  depths  (ft)  at  the  cell  corners 
of  the  boundary-fitted  grid  from  FILE  15  (ITRAN=2).  The 
coordinates  are  then  multiplied  by  the  scale  factor,  XMAP, 
and  divided  by  XREF  to  make  them  nondimensional. 
Subroutine  BJINTR  is  called  to  provide  the  coordinate 
derivatives  needed  to  compute  the  metrics  of  the 
transformation. 


BJINTR  Computes  various  coordinate  derivatives  and  sets  the  water 
depths  HU(I,J)  and  HV(I,J)  on  the  faces  of  each 
computational  cell. 
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CH3DIH  Prints  water  depths,  if  requested  by  input  data.  Also,  the 
water  depths  are  made  nondimensional  by  dividing  by 
ZREF. 

CH3DV 1  Sets  up  flags  for  different  grid  cells  to  indicate  land  and 

water  for  later  use  in  CH3DV2  for  flow  visualization. 

CH3DV3  Similar  to  CH3DV1  for  later  use  in  CH3DV4  for  salinity. 

CH3DV5  Similar  to  CH3DV1  for  later  use  in  CH3DV6  for 

temperature. 

CH3DND  Normalizes  several  variables  and  parameters,  such  as  the 
Eckman  number,  Rossby  number,  time-step,  etc. 

CH3DII  Sets  up  the  arrays  of  boundary  flags  that  indicate  the  nature 
of  computational  cell  boundaries.  In  addition,  arrays 
controlling  the  computation  of  the  convective  terms  in  the 
momentum  equations  and  the  water  surface  cross-derivative 
terms  are  set  up.  One-dimensional  channel  cells  are 
identified. 

CH3DIF  Initializes  various  variables  for  a  cold  start  run  and  opens 

time  series  output  files  for  elevation,  velocities,  salinity,  etc. 
as  well  as  print  and  snapshot  files.  The  hot  start  capability 
is  not  operational. 

CH3DIV  The  arrays  related  to  water  surface  cross-derivatives  created 
in  CH3DII,  contain  logical  values.  Those  arrays  are  used  in 
this  subroutine  to  create  arrays  containing  numerical  values. 
These  arrays,  i.e.  AFV1(I,  J),  etc.  are  used  to  control 
computation  of  not  only  water  surface  cross-derivatives  but 
other  variables  as  well. 

CH3DWS  Controls  the  reading  of  either  wind  speed  or  wind  stress.  If 
the  wind  speed  is  read,  the  stress  is  computed  from  Garratt’s 
equation.  ENTRY  CH3DWT  controls  the  reading  of  time- 
varying  values  and  computations. 

The  subroutines  above  are  called  from  CH3D  in  the  sequence  given. 

Before  calling  CH3DM2,  which  controls  the  computations,  the  initial  salinity 
field  is  read  from  FILE  74.  The  initial  temperature  field  is  read  from  FILE 
17  and  made  dimensionless. 

CH3DM2  Final  subroutine  called  from  CH3D.  All  subroutines 

controlling  the  actual  3D  computations  are  called  from  this 
subroutine  in  the  order  they  appear  below. 
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CH3DDP  Computes  the  total  water  depths  from  the  latest  water  surface 
elevation  field.  ENTRY  CH3DDM  sets  total  water  depths  at 
the  intermediate  time  level  M  and  ENTRY  CH3DDN  sets 
total  water  depths  at  time  level  N. 

CH3DTK  Reads  equilibrium  temperatures  and  surface  heat  exchange 
coefficients  from  FILE  19  and  then  transforms  them  into 
nondimens ional  form.  ENTRY  CH3DTB  controls  the 
reading  of  time-varying  values  and  interpolation. 

CH3DRJ  Reads  river  inflows  from  FILE  13.  ENTRY  CH3DRV 

controls  the  reading  of  time-varying  values  and  interpolation. 

CH3DTI  Reads  and  initializes  tidal  boundary  conditions  from  FILE 
16.  ENTRY  CH3DTD  updates  boundary  values. 

CH3DSAI  Reads  salinities  and  temperatures  at  tidal  boundaries  from 

FILE  76.  ENTRY  CH3DSAV  controls  the  reading  of  time- 
varying  values,  interpolation,  and  conversion  to 
nondimensional  form. 

CH3DTEI  Reads  temperatures  at  river  inflow  boundaries  from  FILE 
78.  ENTRY  CH3DTEV  controls  the  reading  of  time- 
varying  temperatures  and  interpolation. 

WQMOUT  This  is  the  main  water  quality  module.  It  reads  information 
on  cell  corner  points  and  water  quality  boxes  and  computes 
box  lengths,  volumes,  and  flow  face  areas.  ENTRY 
WQMCVOL  computes  the  volumes  for  the  first  time  and 
ENTRY  WQTVD  computes  time-varying  information  and 
stores  it  for  later  input  to  the  Water  Quality  Model. 

PROCl  This  subroutine  called  by  WQMOUT  performs  Stokes’  drift 
computations. 

PTOUT  This  is  the  main  particle  tracking  module.  It  saves  time- 
invariant  information  on  grid  characteristics  first.  ENTRY 
PTCVOL  designates  the  output  file  number  and  ENTRY 
PTTVD  saves  time-varying  information  on  wind  stress, 
water  depth,  surface  velocity,  horizontal  velocity  and 
acceleration,  and  vertical  velocity  for  later  use  by  the 
Particle  Tracking  Model. 

EK2,  EK2V  These  subroutines  called  by  PTOUT  compute  the  Ekman 
drift. 

At  this  point,  the  loop  over  time  is  entered  in  CH3DM2  and  each  of  the 
subroutines  below  is  called,  usually  each  time-step. 
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WQMCVOL  ENTRY  in  WQMOUT  for  computing  initial  box  volumes. 


PTCVOL 

CH3DDE 

CH3DED 

CH3DWT 

CH3DTB 

CH3DRV 

CH3DTEV 

CH3DTD 

CH3DSAV 

CH3DDN 

CH2DXY 

CH3DDP 

CH3DXYZ 


CH3DDI 


CH3DSA 

CH3DTE 


ENTRY  in  PTOUT  for  designating  the  output  file. 

Computes  water  densities  using  Equation  8.  The  baroclinic 

terms  in  the  momentum  equations  are  then  evaluated. 

Sets  up  arrays  that  are  then  used  in  subroutine  CH3DEZ  for 
the  computation  of  eddy  viscosity  and  eddy  diffusivity 
coefficients. 

ENTRY  in  CH3DWS  for  reading  time-varying  wind  data. 

ENTRY  in  CH3DTK  for  reading  time-varying  equilibrium 
temperature  and  heat-exchange  coefficient. 

ENTRY  in  CH3DRI  for  reading  time-varying  river  flows. 

ENTRY  in  CH3DTEI  for  reading  time-varying  temperatures 
at  river  inflow  boundaries. 

ENTRY  in  CH3DTI  for  reading  time-varying  tide  data. 

ENTRY  in  CH3DSAI  for  reading  time-varying  salinity  and 
temperature  profiles  at  tidal  boundaries. 

ENTRY  in  CH3DDP  for  assigning  total  water  depths  at  time 
level  N. 

Computes  the  vertically  averaged  flow  field  from  the 
vertically  averaged  equations  of  motion. 

Computes  the  total  water  depths  at  time  level  N-bl,  using 
the  water  surface  field  computed  in  CH2DXY. 

Computes  the  3D  velocity  field.  Mass  conservation  is 
ensured  by  forcing  the  vertical  sum  of  the  horizontal 
components  of  the  3D  velocity  to  match  the  vertically 
integrated  values  computed  in  CH2DXY. 

Computes  the  convective  and  diffusion  terms  in  the 
momentum  equations  using  the  most  recent  computation 
results  from  CH3DDP  and  CH3DXYZ.  These  terms  are 
then  employed  at  the  next  time  step  in  CH2DXY  and 
CH3DXYZ. 

Computes  the  salinity  field. 

Computes  the  temperature  field. 
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CH3DBL  Checks  the  water  surface  elevations  for  the  program 
"blowing  up". 

WQTVD  ENTRY  in  WQMOUT. 

PTTVD  ENTRY  in  PTOUT. 

CH3DOT  Controls  the  output  printed  and/or  written  to  files  for 

plotting.  Output  is  in  terms  of  physical  dimensional 
variables.  Subroutine  CH3DC1  is  called  with  ENTRIES 
CH3DC2,  CH3DC3,  CH3DC4,  CH3DC5,  CH3DC6, 
CH3DC7,  CH3DC8,  CH3DC9,  CH3DCA,  CH3DCC, 
CH3DCD,  and  CH3DCE.  Each  is  described  below. 

CH3DC1  Provides  dimensional  water  surface  elevations. 

CH3DC2  Provides  dimensional  physical  vertically  averaged  velocity  in 
jc-direction. 

CH3DC3  Provides  dimensional  physical  vertically  averaged  velocity  in 

y-direction. 

CH3DC4  Provides  dimensional  physical  horizontal  velocity  component 
in  ^-direction. 

CH3DC5  Provides  dimensional  physical  horizontal  velocity  component 
in  y-direction. 

CH3DC6  Provides  dimensional  physical  vertical  component  of  3D 
velocity. 

CH3DC7  Provides  salinity. 

CH3DC8  Provides  dimensional  temperature. 

CH3DC9  Provides  dimensional  physical  magnitude  and  direction  of 

horizontal  velocity. 

CH3DCA  Provides  dimensional  physical  horizontal  components  of  3D 
velocity  at  the  centers  of  cells. 

CH3DCC  Provides  dimensional  water  density. 

CH3DCD  Provides  dimensional  vertical  eddy  viscosity. 

CH3DCE  Provides  dimensional  vertical  eddy  diffusivity. 
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After  CH3D0T  the  following  subroutines  for  visualization  are  called  in 
CH3DM2  in  the  order  shown: 


CH3DV2  Saves  time-varying  flow  information. 

CH3DV4  Saves  time-varying  salinity  information. 

CH3DV6  Saves  time-varying  temperature  information. 

In  subroutine  CH3DOT,  the  following  files  are  created  for  use  in 

generating  time  series  plots,  vector  plots,  or  contour  plots. 

FILE  21  For  time  series  plots  of  dimensional  water  surface  elevation 
at  specified  horizontal  locations. 

FILE  22  For  time  series  plots  of  dimensional,  Cartesian  horizontal 
velocities  (x  and  y  directions)  at  cell  centers  at  specified 
horizontal  locations. 

FILE  23  Geometry  of  study  area  (needed  for  plotting  snapshots 
or  contours). 

FILE  24  For  velocity  vector  plots  and  contour  plots  of  surface 
elevation,  salinity,  temperature,  etc. 

FILE  25  For  time  series  plots  of  discharges  at  specified  horizontal 
ranges. 

FILE  31  For  time  series  plots  of  salinity  at  specified  horizontal 
locations  in  all  layers. 

FILE  34  For  time  series  plots  of  temperature  at  specified  horizontal 
locations  in  all  layers. 

FILE  35  For  time  series  plots  of  vertical  eddy  viscosity  at  specified 
horizontal  locations  in  all  layers. 

FILE  36  For  time  series  plots  of  vertical  eddy  diffusivity  at  specified 
horizontal  locations  in  all  layers. 

FILE  37  For  time  series  plots  of  density  at  specified  horizontal 
locations  in  all  layers. 


The  model  creates  the  following  additional  output  files. 
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FILE  42 
FILE  43 
FILE  44 
FILE  70 
FILE  71 
FILE  96 


For  visualization  of  flow. 

For  visualization  of  salinity. 

For  visualization  of  temperature. 

Geometric  data  and  initial  results  for  particle  tracking. 
Time-varying  data  for  input  to  the  Particle  Tracking  Model. 
Time-varying  data  for  input  to  the  Water  Quality  Model. 
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4  Demonstration  of  the  Setup 

of  Input  Files 


To  demonstrate  the  setup  of  input  files  (see  the  appendixes),  portions  of  the 
input  files  are  presented  for  an  application  in  which  the  hydrodynamics  of  the 
New  York  Bight  during  April  1976  was  simulated.  This  application 
represents  the  beginning  of  the  long-term  (6-month)  model  validation.  Results 
from  this  application  are  presented  by  Schefftier  et  al.  (1994). 


INCLUDE  Files 


As  previously  indicated,  two  INCLUDE  files,  nyb.inc  and  chSd.inc,  are 
needed  for  compiling  the  model.  Of  these,  nyb.inc  is  used  to  define  variable 
array  dimensions  in  the  model  for  running  a  particular  application,  e.g..  New 
York  Bight  6-month  simulation.  It  may  be  necessary  for  the  user  to  change 
some  of  the  parameters  in  nyb.inc.  It  is  listed  in  Appendix  A  (page  A-1). 
Parameters,  whose  values  are  set  in  nyb.inc,  are  used  to  dimension  arrays  in 
COMMON  blocks  in  ch3d.inc  and  other  arrays  in  the  model.  All  the 
parameters  appearing  in  nyb.inc  are  defined  below.  Of  these,  ICELLS, 
JCELLS,  UMAX,  and  KM  have  to  be  set  exactly.  The  others  can  be  greater 
than  or  equal  to  what  is  needed. 


ICELLS 

JCELLS 

UMAX 

KM 

NSTATS 

NTIDES  =  11 
NRIVRS 
NCNST  =  37 


NBNDS 

NBARRS 


:  Number  of  grid  cells  in  the  ^-direction 
:  Number  of  grid  cells  in  the  77-direction 
:  The  greater  of  ICELLS  and  JCELLS,  plus  1 
;  Number  of  u-layers  in  the  vertical 
:  Maximum  number  of  gauge  stations  where 
information  will  be  saved 
:  Not  used 

:  Number  of  river  boundaries  used 
:  Maximum  number  of  tidal  constituents  used 
(set  to  37)  -  used  only  if  tidal  signals 
were  generated  using  constituents  -  not 
operational. 

;  Number  of  open  water  boundaries  used 
:  Number  of  interior  thin-wall  barriers 
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NPRWIN 

Number  of  print  windows  for  printing 
model  results 

NSNAPS 

Number  of  snapshot  windows  where 
information  is  saved 

NRANGS 

Number  of  ranges  where  discharge 
information  is  saved 

NTIDFN 

Number  of  tide  functions  used 

NTIDBN 

Number  of  tidal  boundaries  used 

NTIDPT 

Maximum  number  of  values  in  the  input 
tide  functions 

NWINDS 

Number  of  input  time-varying  winds  used  for 
interpolating  winds  over  grid  (special  for 

New  York  Bight  application) 

NROWS 

Maximum  number  of  computational  chains  used 
in  ^-direction 

NCOLS 

Maximum  number  of  computational  chains  used 
in  ij-direction 

KROWS 

Larger  of  NROWS  and  NCOLS,  plus  1 

NX8PTS 

Number  of  one-cell  wide  channel  cells  in 
^-direction 

NY8PTS 

Number  of  one-cell  wide  channel  cells  in 
jj-direction 

ISSMAX  =  10 

Not  used 

NPLPTS  =  75 

Not  used 

NSDMAX  =  10 

Not  used 

SPVAL 

A  small  value  to  which  the  vertical  eddy 
coefficients,  etc.  are  set  as  a  default 

In  addition,  the  following  parameters  are  used  in  connection  with  the  water 
quality  module  in  CH3D.  They  should  be  changed  only  if  the  Water  Quality 
Model  grid  is  changed. 

NBP  :  Number  of  water  quality  boxes  in  plan  view. 

NFP  :  Number  of  flow  faces  in  plan  view. 

A  portion  of  the  file  chSd.inc  is  listed  in  Appendix  A  (page  A-2).  This  file 
does  not  have  to  be  changed  from  application  to  application,  but  remains  the 
same. 


Basic  Control  Data 

Appendix  B  lists  the  input  variables  and  their  format  for  the  primary  input 
file  FILE  4.  Appendix  C  gives  a  list  of  all  the  input  files.  Input  data  in  FILE 
4  for  the  application  mentioned  are  listed  in  Appendix  D.  These  data  are  in 
response  to  read  statements  presented  in  Appendix  B.  As  can  be  seen,  the 
computational  time-step  (DT)  is  150  sec  with  a  total  of  17,280  time-steps 
(30  days)  simulated  (IT2=  17280).  Both  temperature  and  salinity 
computations  are  made  aTEMP=  -2  =ISALT),  with  the  initially  specified 
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fields  frozen  for  the  first  2,880  time-steps  (ITSALT=2880).  A  time-varying 
wind  (IWIND=5)  is  specified.  There  is  one  river  (NRIVER=1),  the  Hudson. 
Eight  tidal  boundaries  (TIDBND  =  8)  and  eleven  tide  functions  (TIDFNO=ll) 
are  used  to  describe  the  time-varying  tides  at  the  boundaries  by  linear 
interpolation. 


Freshwater  Inflows 


Daily-averaged  freshwater  flows  specified  at  the  Hudson  River  boundary 
are  read  from  FILE  13.  They  are  given  in  Appendix  E  for  the  first  few  days 
of  the  simulation.  Input  fields  are  day,  hour,  location  of  the  river  boundary, 
and  the  discharge  (cfs). 


Wind  Speed 


Wind  data  are  read  from  FILE  14.  Since  IWIND  =  5,  time-varying  wind 
velocity  components  (m/sec)  are  read  and  wind  stresses  are  computed  using 
Equation  11.  For  the  April  1976  simulation,  three  wind  fields  are  read 
(NWINDS  =  3).  These  correspond  to  hindcast  data  obtained  from  three 
offshore  stations  of  WES  Wave  Information  Study.  Inside  CH3D-WES,  wind 
velocity  components  for  each  computational  cell  are  assigned  by  linear 
interpolation  from  the  values  read  in  for  the  three  offshore  stations.  As  can 
be  seen  in  Appendix  F,  each  line  of  input  in  FILE  14  consists  of  the  day, 
hour,  and  the  x  and  y  components  of  the  wind  velocity  for  each  of  the  three 
wind  stations.  Wind  data  at  the  three  stations  are  presented  in  Appendix  F 
only  for  the  first  few  days  in  April  1976.  Wind  stresses  computed  in  the  x 
and  y  directions  are  transformed  in  the  model  into  contravariant  components. 


Grid  Coordinates  and  Water  Depths 


The  X-  and  y-  coordinates  of  the  cell  corners  of  the  computational  grid 
(Figure  2)  and  the  corresponding  depths  (in  feet)  are  read  from  FILE  15.  The 
coordinates  were  generated  using  the  grid  generation  program  EAGLE 
(Thompson  1987a,  1987b).  Landlocked  points  were  assigned  coordinate 
values  of  9.0  x  10  Cell  water  depths  were  obtained  by  interpolation  from 
bathymetric  data  digitized  from  National  Ocean  Service  nautical  charts  for  the 
New  York  Bight  region.  The  format  of  these  data  for  a  few  points  is 
illustrated  in  Appendix  G.  The  first  line  shows  the  title  and  the  second  line 
the  number  of  points  in  the  ^  and  rj  directions.  The  x-  and  y-  coordinates  and 
depth  of  each  point  are  then  input  starting  with  row  1  and  moving  from  left  to 
right,  one  point  per  line. 
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Tabular  Tides 


The  tabular  tide  data  are  read  from  FILE  16.  There  are  8  tidal  boundaries 
and  the  tide  is  interpolated  at  the  boundary  cells  using  1 1  tide  functions,  as 
defined  in  the  basic  input  file.  Each  string  of  input  contains  the  month,  day, 
year,  hour  and  minute  along  with  the  corresponding  values  of  the  water 
surface  elevation  (cm)  at  the  11  tidal  function  point  locations.  Appendix  H 
illustrates  the  form  of  these  data  for  the  first  few  hours  of  April  1976. 


Initial  Temperature  and  Salinity 

FILES  17  and  74  are  identical  in  form.  FILE  17  contains  the  initial 
temperature  field,  whereas  FILE  74  contains  the  initial  salinity  field.  These 
files  were  generated  by  a  computer  program  which  uses  a  few  observed  values 
to  assign  values  to  the  individual  cells  in  each  of  the  sigma  layers.  The 
resulting  field  is  then  smoothed  in  each  horizontal  direction  in  each  layer  and 
the  smoothed  field  is  written  to  either  FILE  17  or  FILE  74.  Appendix  I 
presents  the  initial  temperature  field  for  sigma  layer  1 . 


Surface  Heat  Exchange  Information 


Daily  averaged  equilibrium  temperatures  and  surface  heat  exchange 
coefficients  computed  from  meteorological  data  at  the  John  F.  Kennedy 
Airport  are  read  from  FILE  19.  As  indicated  in  Appendix  J,  each  input  data 
line  consists  of  the  day  and  hour  since  the  initiation  of  the  simulation,  the 
equilibrium  temperature  (deg  C),  and  the  surface  heat  exchange  coefficient 
(cm/sec). 


Tidal  Boundary  Salinity  and  Temperature 


At  tidal  boundaries,  the  time-varying  vertical  distributions  of  salinity  and 
temperature  are  read  from  FILE  76.  As  illustrated  in  Appendix  K,  these  data 
are  grouped  in  the  following  manner.  The  day  and  hour  since  the  initiation  of 
simulation  is  input  on  one  line.  It  is  followed  by  a  line  giving  the  (I,J) 
location  of  the  tidal  boundary  cell  and  the  corresponding  vertical  salinity 
distribution.  The  next  line  gives  the  same  (I,J)  location  and  the  corresponding 
vertical  temperature  distribution.  These  distributions  consist  of  values  in  each 
sigma  layer,  starting  with  the  surface  layer  and  going  down  to  the  bottom. 
This  input  is  repeated  for  all  the  tidal  boundary  cells.  Appendix  K  gives  the 
input  for  the  61st  day  of  simulation  for  a  few  boundary  segments.  Note 
information  has  to  be  repeated  for  cells  common  to  two  adjacent  boundary 
segments. 
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River  Temperature 

River  temperature  data  are  read  from  FILE  78  (APPENDIX  L).  In  the 
present  case,  there  is  only  one  river  boundary,  the  Hudson.  The  values  are 
estimated  from  available  U.S.  Geological  Survey  gauging  station  data.  The 
first  line  contains  the  day  and  hour  since  initiation  of  the  simulation.  Next  for 

each  of  the  NRIVER  rivers  indicated  in  the  input,  the  I  and  J  cell  indices  are 
given  followed  by  the  temperature  in  different  layers,  starting  from  the  surface 
and  going  down  to  the  bottom.  This  is  done  in  the  order  URSTR  to  IJREND 
for  each  river.  The  same  procedure  is  followed  for  successive  times. 
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5  Summary 


The  primary  purpose  of  this  report  is  to  serve  as  a  user  guide  for  the  New 
York  Bight  3D  hydrodynamic  model.  After  a  brief  introduction  to  the  New 
York  Bight  study  region,  the  New  York  Bight  feasibility  study  and  the 
computational  grid  used,  the  report  describes  the  main  features  of  the  CH3D- 
WES  hydrodynamic  model  selected  for  the  study.  In  Chapter  2,  the  basic 
governing  equations  are  given  followed  by  the  boundary  and  initial  conditions 
employed. 

The  report  outlines  the  structure  of  the  computer  model  in  Chapter  3, 
listing  the  names  of  the  various  subroutines,  their  functions,  and  the  calling 
sequence.  This  should  help  users  who  are  interested  in  following  the  logic  of 
the  model.  Also  listed  are  various  output  files  created  by  the  model  and  their 
contents. 

Chapter  4  of  the  report  describes  the  INCLUDE  files  needed  for  defining 
grid  dimensions,  COMMON  blocks  ,  etc.,  and  discusses  the  basic  control 
data,  and  various  input  files  for  a  specific  application  for  the  New  York  Bight 
study  region  covering  the  month  of  April  1976.  A  partial  listing  of  the 
INCLUDE  files  is  provided  in  Appendix  A.  The  format  of  the  primary  input 
file  is  given  in  Appendix  B  followed  in  Appendix  C  by  a  list  of  all  the  input 
files  required.  The  user  guide  function  of  the  report  is  enhanced  by 
demonstrating  the  various  input  files  in  Appendixes  D-L  for  the  specific 
application. 
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Appendix  A 

Listing  of  INCLUDE  Files 


c 

C-  ■  INCLUDE  FILE  FOR  NY  BIGHT  19935  revised  by  Rao  7/29/93 

C 

C 


PARAMETER  <  ICELLS  =  76. 

XELLS  =  45, 

IJfttX  =  77  ) 

PARAMETER  (  KM  =  10, 

KIM  =  KM  -  1 

) 

PARAMETER  (  NSTATS  =  40. 

KTIDES  =  11, 

MtIVRS  =  9  ) 

PARAfETER  (  NCNST  =  37, 

NBNDS  =11, 

N»«RS  =  9  ) 

PARAMETER  <  lERHIN  *  20, 

fSWS  =  20, 

f«AN6S*20  ) 

PARAfETER  (  KHOFN  =  12. 

KTIMN  =  11, 

NTIDPT  =  10000  ) 

PARAfETER  <  fWINDS  =  3 

) 

PARAfETER  (  fTOHS  =  120. 

NCOLS  =  120, 

KROHS  =  121) 

PARtfETER  (  NX8PTS  =  10.  NY8PTS  *  10  ) 

C 

PARftfETER  (  ICl  =  ICflLS  ♦  1.  XI  =  XEliS  ♦  1  ) 

PflRflfETER  <  in  =  ICELLS  +  2.  JH  =  XELLS  ♦  2  ) 

PARAMETER  (  ISSlttl  =  10.  NPIHS  »  75  ) 

PARAMETER  {  NSDMAX  =  10.  SPW.  =  .1234E-6  ) 

C 

added  PARAMETERS  frt*  Hater  Buality  Module  !  Ru  7/29/93 
C***  Caution  1  Do  not  change  inless  Hater  Quality  Model  grid  changes  !!! 
PARAfETER  (NBP  =  3000.  ffP  «  6000) 

C*«*«« 

c 
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A1 


c«***  cti3d.inc  for  KY  Bight  project;  last  revised  by  Rao  7/29/93  -  includes 
c  particle  tracking  variables;  added  flags  for  saving  info  for  water  quality 
c  (IHQ)i  particle  tracking  (IPT)i  and  visualization  (IVIS). 

C — INOUIie  FILE  =  CH3D.  INC  Last  revised  on  2  fttY  91  Klfl 

C 

COttlON  /BLKOl/  HS{0:IH.0:JH1,  HU«0:Itt.0:J«),  HV(0;IM,0:J!1)i 

1  HSD<o:itt.o:Jn) 

comoN  /HJ(02/  ui(o:i«,o:»).  uih(o:ih,o:jh),  uin(o:ik,o:j«) 
catON  /BLKCO/  vi(0:i«.o:j«).  vin(0:iM,o:3«),  vin(o:ik,o:») 
coftiQN  /BU04/  s(o;i«.o:3«i.  ai(o:iH.o:J«).  SN(0:i«,o:») 

C 

ctrnw  /BLKOS/  AHUU(0:i«.o:ai).  fiHJH(o:i«.o:J«)i  ahlk(o:i«.o:j«). 

1  fWV(o:iH.o:«)i  AiMi(o:iH.o:jn)t  mvn(0:ih.o:jn). 

2  WfiS(o:lK.o:»).  fiHSN(0:lH.o:J«).  tfS«{0:lK.0:»). 

3  AHUI(0:IH.0:3K}.  fWI(0:IH.0:JH>. 

4  AKSU(o:iti.o:m>Kn)>  <^(o:in.o:ffi.Kn> 

C 

ctmw  /BLKo;/  tfVi<o:i«,o:»).  afv2{o:ih.o:W).  tfV3(o:i«.o:JH). 

1  tfV4(o:i«.o:w).  tfui(o:i«.o:»).  tfU2(0:iii.o:3«). 

2  AFU3(o:iH.o:}n).  AFm(0:in.o:3ii) 

c 

amw  /Buoe/  EEXdci.Kii.  dosycici.jci), 

1  Diandci.xi),  mDYdci.xi). 

2  DHSXdCl.JCl),  DHSYdCl.Xl). 

3  DHTXdCl.Xl).  DHTYdCl.Xl). 

4  DHUXdCl.XI).  DHUYdCl.JCl) 

C 

COttON  /BU09/  OllXdCl.Xl).  DlIYdCl.Xl) 

comm  /BUdo/  xoidci.xn.  XD2dci.xi)>  XD3(ici.xi)> 

1  YDldCl.XI).  YD2dCl.Xl).  YD3dCl.Xl) 

comoN  /BLKii/  RiRHOiOdn.oiJn.KH).  DiRHQ(o:in.o:jn.XH} 

COmON  /HJ(12/  THETA 
comm  /ADD77/  KU.  KV 

comm  /AD078/  BffiU(o:in.o:JH.XH).  Bffiv(o:in.o:m.KH> 


comm  /BLK16/  XUdCll.  XSdCll.  YV(Xl).  YS(Xl).  ALXREF,  ALYREF 
comm /suci?/ NS(o:iK.o:m}.  irohinrous).  iukmioks). 

$  IU2(NRO»5).  lEMOflOUS).  KlOH 

comm  /BUdB/  iis(o:iH,o:m).  xol(ncols).  jvknccls), 

%  JV2(NC0LS),  jmoeas).  ncol 

comm  krokckrousi.  xukxrous),  kuzikrous).  xkol 

COMNm  /BLK19/  NR(0:I«.0:JH).  KR(0:IN.0:JH) 
comm  /BU20/  DTi  DIT,  DX.  DY,  DTDX,  DTBY.  DXI.  DYI,  D2XI.  D2YI. 
i  DELT 

comm  /BU22/  Dox(o:i«.o;m),  DeDY(o:i«.o;m). 

1  fiix<o:ib.o:jh).  fiiy(o:in.o:j?1). 

2  FH)X(o:iK.o:m).  FH3Y(Q:in.o:jn}. 

3  coRX(o:iK.o:m).  CQRY(o:in.o:m). 

s  xxaR(o;in.o:m.xH)t  YKORioiin.oijniXHit 

4  TBX(0:IK.0:3H}.  TBYiOtlHtOlJm.  WO 

comm  fBJW  imDIR(fRIVRS).  IJRRONCNRMB),  IJRSTRORIWS). 

1  IXENDitSUVRS).  ORIVER<WIWS.IJnAX)i  WI\ER. 

2  IDAYA>  IHOURAt  IDAYB7  IKXiS 

comm  /BU(24/  IJKJIR(NBARRS)t  IJBRONClfflflRRSli  IJBSTROfflARRSlt 
$  IJBEKD(NBARRS)t  NSAR 

comm  /BU(25/  IJDIR(HBNDS>t  IJROW(NBM)S)i  IJSTRT(HBNDS)t 
$  IJEIQKNBNDS).  HIDE.  IXINE.  NCG.NCONST 


A2 
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c 

comoN  mmi  amphcxellSiKtides).  am^(jcells<ktii£S). 

1  Wf>S(ICaLS,KTIIES).  AMPNdCELLS.KTIIlES). 

2  PrtKmLS.KTIDEB).  PHE<X£LLS,KTIDeS), 

3  PKS(ICEliS.MTIDEB).  PHNdCEULS.MTIDES). 

4  CrtdJCELLS.KTIDES),  Cfl£(XEliS.NTIDES), 

3  CftSdCELLS.KTIDeS).  CfiNdCEllS.NTIDES). 

6  P(MTIDES) 

C 

OmJN  /BU27/  CRIVRA(NRIVRS.IJHAX).  6RIVRB(NRIVRS,IJKa) 
COmON  /HJ(28/  IST(feTATS).  JSTWETATS).  NSTA,  tFRES 

C - 

OWACTERMa  STATID.  STATSA,  STATS 

C - 

COmON  /BLI(29/  STATID(NSTATS).  STATSA (NSTATS),  STATS(N5TATS> 


COrtCM  /BLX30/  ITBRTdO).  RSPACdO).  ISPACdO).  JSPACdO). 
t  imXl.  ITBRKI 

COmON  /BU31/  ITRAN.  ISO 
comoN  /ax22/  TYCOtiH.o:^). 

1  TX1(0:I«.0:J«),  TY1{0:IM.O:J1). 

2  'nc2(o:iti.o:jH).  TY2{o:ih.o:Jti)> 

3  TAUX(onH.o:»).  TAuy(o:iii,o:J?i). 

4  TALai(o:i«.o:j«).  taux2{o:i«.o:j«). 

5  tauyi(0:im.o:j«).  tauy2(o:im,o:jk), 

6  MNDXlMtlNDS).  »{0X2(NUINDS) . 

7  WffiYlUMIKDS).  i9CY2<NH!N!S).t(SPD(ICl.Xl)i 

8  lOAYlt  IHOURliTIKEli  IDAY2i  IHQUR2i  TII£2. 

9  RHOAIR,  ICtlNVF.TXF(0:lW,0:a).7YF(0;I«,0:Jn), 

.  txfl<o:im.o:»).tyfl(o:im.o:jh) 

ccmw  /BLJC33/  ITESTf  IPA.  IPS.  ID.  JPA.  Jra.  JD 
CmCN  /BLK34/  RWldCELLS.XELLS).  COR.  GR.  CTB.  ROI.  CZNJLT 
OmiN  /B(J(3S/  I6I.  IGH.  I6T.  IGS.  IGU.  IGH.  IGC.  IGQ.  IGP 
COnraN  /BUC36/  IFI.  IFD.  IFD5. 

1  IW.  IPT.  IVIS 

COtHM  /RXJ7/  IT.  IDAYO.  IHOURO.  ItllNO.  ISECO.  HI.  IT2.  TIHE. 

1  Tiro.  nro.  itws.  itsalt.  Titoc.TirEA.  tiio. 

2  Tiro.  Tiro.  IDAYT1.IH?T1.IDAYT2,IHRT2. 

3  TItETl.TirET2.  TEP2.  TEK2.  ITPTS 
OmM  /EU(38/  iNiro 

OnCN  /EU(37/  IBTrt.  HI.  K2.  IKID.  HHIN.  SSSO.  UttP 
aXHM  /ELK39A/  AldJtWO.  A2d^).  A3dJtW(}.  MdJtVUC). 

$  BldJTUX).  B2dma}.  BSdJHAX).  84d3l1AX}. 

i  M(imi),  BSd^).  ccdJKAX).  indJnAX) 

c - 

CHARACTERS  NBXdCELLS.XELLS).  KBYdCEUS.XflLS) 

C - 


COrtfM  /BUC40/  NBI.  NBY 
C 

COtOi  /PFUe/  IPSU.  IPUIH.  IPVIH,  IPSSM.  IPUU.  IPVN. 

1  IPSS.  IPUIG.  IPVIG.  IPSSG.  IPU6.  IPV6. 

2  IPW.  IPSAM.  IPTEH.  IPROH. 

3  1PM6.  IPSA6.  IPTEB.  IPR06 

C - 

COttfM  /TRN017  6D(0:in.0:W.3) 

cottiM  mm/  6U(0:in.o:3M.3).  612(o:ih.o:jh.3). 

$  G22(0:in.o:JH.3) 
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COmON  /TRN03/  Hll(0:lMiO:J«,3).  H12(0:I«.0:J«.3). 

♦  H22{0;i«.o:j«.3) 

amjN  mm/  Diii(o:iN,o;ai.3).  Dii2(o:iN.o:jni3). 

<  D122(0:IM,0:».3) 

Ctmn  /7W»5/  D211<0:IH.0:».3),  D212(0:IM.0:J«,3). 

♦  I>222(0:IM.0;W,3) 

ctncN  mm//  fxii(0:im.o:jh),  n(i2{o;i«.o:JH). 

♦  FX2i(0:iH,0:J«).  FX22(0:l«,0:J«) 

coffON  mme/  fyii(o:im.o:ji),  fyi2{o:i«.o;j«), 

«  FY21(0:iH,0:»).  FY22«0:l«,0:W) 

COffON  /TRN09/  XlU(ICl.Xl)iX2U(ICl,JCl),yiU(IClf JCI>iY2U(ICliJCl) 
COMMON  /TFNIO/  X1V(IC1.X1).X2V(IC1,X1),Y1V(IC1,X1).Y2V(IC1.X1) 
COMMON  /TRNll/  XlC(ICliXl)iX2C(ICl,Xl)iYlC(ICliXI)ty2C(ICl,Xl) 
COMMON  /TKN12/  XI  {IC1.X1).X2  (ICI,X1),Y1  (IC1.XI).Y2  (ICl.XI) 

COmON  /flfiEAOl/  XadCl.Xll.  XCXXdCl.XDi  XCYYdCliXl). 

♦  XCXXYYdCl.Xl) 

COmON  /fiREA02/  YCTdCl.Xll.  YCXXdCl.Xl).  YCYYdCl.Xl}. 

<  YCXXYYdCl.Xl) 


COmON  /fflfflOl/  FMUdCl).  FMV(Xl),  FMSdCl).  FKSV(Xl) 

omw  /fmw  S6(XM).  dzz(xm).  z(xm).  dza(x«),  dzb(xm) 
coftw /AWXKJ/  u{o:in.o:3M.iMi.  UN(o:iM.o:Jn.xw,  un(0:i«i0:jm.km) 
coftON  /flD004/  v(o;iM.o:».xn).  vn(o:im.o:jii.xm),  vh(o:ih.o;jm.kn) 

COMMON  /ADD05/  N(0;IH.O:JM.XM).  WdO;IH.O:ai,XK),  NN(0:iM.0;JMtXM) 
1  .NT(o;iM.o:JM),u3AV<o:iM.o:JM,x«).V3W(o;iM.o:JM,Kn) 

comoN  /ADD06/  FBC(o:iM.o: Ji)t  fbcv(o;im.o:ji). 

4  FI  (o:im,o:ji,km>.  fj  «o:im,o:jm,k«). 

1  Fra(o:iHiO:»,KH),  FDABC(o:in.o:3nii!M), 

2  FRQ(0:IH.0:»,K«).  F8ABC(0:i«,0:;»1,KH). 

3  FXYZ(0:IM,0:3M,KM) 

COtfCN  /ADD07/  C  WllHiOlJM.XH),  CAl0;Ih)0: J1),  CEQ(0:i(1.0:3n). 

1  T  (0:IM,0:3M,K«),  TN  W:  111,0;  JI.XM). 

2  Sfl(o:in.o:»,KM).  sw(o;im.o;3M.km), 

3  SAa(0:lH.O:JK.KM).XldH),  JS2d«),  FSIH.  SAINT 
COftON  /ADDZA/  IDAYSl,  IHRSl,  TlttSl,  IDAYS2,  IHRS2,  TIMES2, 

1  SAl(NTIDEB,I«AX,XM),  SA2(NnEES,IJMAX,XM), 

2  TEKNTIDES.IWAX.XH),  TE2(NTIDES,IJ1AX,XH) 

COMMON  /m/B/  IDYIEl,  IHRTEl,  nuiH,  IDYTE2,  IHRTE2,  nNTE2, 

1  Ta(f«IVRS.I3MAX.XM),  TE4(M?IVRS,IJKAX,XM), 

2  SA3(NRIVRS,I^.Km.  SA4(M;iV1S,im.Kn) 

COmON  /ADOOB/  R(0:in,0;JM,XH),  RU(0:lM,0:JM,KH),RV(0;lM.0;Jn,KM), 

1  RS(0:iM,o:»),  RtJS(0;iM,o:J«),  rvs(0:im,0:JM) 

COmON  /AMXW/  GA  (0:IM.O:JM.K«),  S  {0;IM.0:JM,iniJ, 

1  GAB(  o:iH,o:JM).  o:iN.o:JM). 

2  RIA(0:IM,0:JM.KH).  RIBIOlIH.OtJM.UU. 

3  VE1JGA(0;IM,0:J1,XH),  VELS{0;IM.0:JM,XM) 

COmON  /AUDIO/  0QQ<0:IH.0:3H.KH}.  5L  (0: 1)1,0: Jlt.Kn) 

comoN  /ADOii/  uix  (o;im,o:3M),  uiy  {o;in,o;jM),  uiC{o;ni,o:JM), 

1  vix  (o:im,o:jm),  viy  (o:iii,o:Jn).  vic(o;im.o:3H), 

2  DC  (o:iM,o:jn,K)i).  vc  (o:di,o:jm,kh) 

COmON  /AD012/  1)RES(0:IH,0:«,KM),  VRESIO:  1)1,0:3)1,101), 

1  NRES(0:IH,0:3M,XH). 

2  uires(o:ih,o:jm),  vires(o;im,o:3)1), 

3  TBXRS(o;ifi,o::m),  tbyrs{o:im,o:ji) 
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comDN  /A0D13/  uiN(o:jn)i  uiE(o:JK).  vis(o:iH)i  vim(o:iH). 

1  SABE(Kn).  SAEM(Km.  SABS(Kn)>  SA&N(Kn). 

2  TBE  (Ktl).  mi  m,  TBS  (KH).  TBN  <Kn)> 

3  CBE  (KH)>  CBW  (KH)>  CSS  (KH).  CBN  (KH) 

COMMON  /ADDH/  TXC  (0:1M.0:JM).  TYC  (0;IM,0:3H), 

1  n2C  ty2C  jo:im,o:3M). 

2  fricx(o:i«.o:jm).  fricy(o:in.o:jm)i 

3  Bzo  <o:in.o:;m).  bzou  (o:in>o:jm) 

c 

COtICN  /INPOl/  XDIST(2*IM+1),  )(SLQPE(2«IH+1), 

1  YDlST(2«Jm-l),  YSLQP£(2*IK+1), 

2  XKftPPA(37,KTII]eS).  )njDN6(NTIIlES) 

COMMON  /IfP02Z  AVI,  AV2.  BZl.  CO,  FXB.  2REFBN,  ZREFTN,  Ot€6A 
comoN  /:teo3/  icci,  icc2.  jcci,  xa,  idi.  id2,  jdi,  jd2, 

1  IGRID,  T80,  ,TSE,  'nH, 

2  TZl,  XDAY,  rJOJR,  XYEW.  XMONTH 

c 

COmON  /AD050/  IVER,  ICON,  lUBO,  !BL.  IBR,  JBM,  JEP.  CREF,  CMAX 

COMMON  /AOOSl/  XREF,  ZREF.  UREF.  K£F.  SREF.  ROO,  ROR.  TO,  TR 

COtON  /AI)052/  lEXP,  lAV.  AVR.  m,  tPA.  AVREF.  AVM,  AVMl 
COmON  /ADD53/  IVia,  IS^T.  ICC,  IFA.  IFB,  IFC 

COmON  /AOOS4/  IRO,  IH,  IMI,  ICI.  IR4.  ItC.ICONC.  INS.  IREAD.  IRON 

COmON  /AD0S5/  DZ,  EH,  EV,  H,  FR,  FR2.  FRD,  FRD2,  RB,  RBV.  BETA, 

1  SZ.  TAUR,  KIS,  IfTU,  inV 

COmON  /ADD56/  IPl,  IP2,  IP3,  IPU,  IPN,  XPA,  XPB,  KD, 

1  IGL,  IGRI,  IGTB,  IRES.  TRES 

COMMON  /ms?/  ITEMP.  TQ(0:IM.0:JM),  BVR.  SI.  S2.  PR.  PRV 
COMMON  /AD058/  NR^GE,  KSTINSTATS) 

COItON  /ADIJ59/  FMl,  FM2,  ZTOP.  BSC.  ZAB,  TRIl.  TRI2.SaJl,  SOC. 

1  Q2R1.  WX 

COMMON  /ADD60/  ISfWl.  ISIE,  ITB 

COmON  /M)D61/  ATS.  ETB.  STB.  VCO.  NK.  ZREFB.  ZRER.  ZOT. 

1  SLHIN,  QQNIN 

comoN  /m62/  ecuT.  icur.  ga^x,  smax,  fzs,  ksmall 

COMMON  /ADD&3/  SSTVOt.  ISTVU.  J9VUC.  UIWX.  lUKU.  JUMAX.  XUmX, 

1  WMAX.  IVHAX.  3VMAX.  XVWX 

COMMON  /AD064/  I)HX(0:IM,0:^).  BKYCOZIN.OIJM) 

COmON  /AOO&S/  IX8(NX6PTS).  lYBINXBPTS).  NX8, 

1  JX8<NY8PTS),  JYBIKYBPTS).  NYB 

COMMON  /ADD6&/  NISS.  ISS(ISSWX),3SS(ISSmX), 

1  NDEPTH  INSDtttX),  TI)EPTH(NSI»«X). 

2  ZDEPTH  (NSOmX.ISEZttX). 

3  ZSAL  (KSDMAX.IS9IAX).  SK(XH.ISSmX) 

comoN  //m,?/  sx(o;im,o:jm).  sy(o:im,o:jmi 

coticN  //me/  acst(ncnst),  pcsTitoen 
cotifM  /AD06Y/  tep(o:ih.o:jh).  TEX(o:iH.o:jn).  HN(o:iN,o:^> 
c 

comoN  misj/  areaioiim.oijm).  dfac(Ici.xi> 
c 

c**f*  added  for  visualiiation  Rao  7/29/73 

INTEGER  VISSTRF,  VISENDF.VISINTF,  VISSTRS,  VISENDS,  VISIKTS. 

1  VISSTRT.  VISENDT.  VISINH 
COmON  /VIS/  VISSTRF,  VISEWF.  VISINTF.  VISSTRS.  VISEWS. 

1  VISINTS,  VISSTRT,  VISEIOT.  VISIMH 

C*M» 
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Appendix  B 

List  of  Input  Data  in  File  4 


DUMMY 

TITLE  Run  descriptor  (Format  A80) 

DUMMY 

IT1,  IT2,  DT,  ISTART,  ITEST,  ITSALT  (218,F8.0,4I8) 

IT1 

Starting  time  step  (  always  set  =  1) 

iT2 

Ending  time  step 

DT 

Computational  time  step  in  sec 

ISTART  =  0 

Cold  start 

>  0 

Hot  start  (not  operational) 

ITEST  =  0 

No  diagnostic  output 

>  0 

Diagnostic  output 

ITSALT 

Number  of  time  steps  after  which  salinity  and 
temperature  computations  are  initiated 

DUMMY 

WPRCRD  (9I8,A8) 

WPRCRD 

Number  of  print  control  lines  which  follow 

DUMMY 

WXCEL1,  WXCEL2, 

WYCEL1,  WYCEL2,  WZCEL1,  WZCEL2,  WPRINT, 

WPRSTR,  WPREND,  WPRVAR  (9I8,A8) 

If  WPRCRD  >  0,  WPRCRD  lines  have  to  be  furnished  below. 

WXCEL1 

Starting  |-cell  index 

XCEL2 

Ending  |-cell  index 

WYCEL1 

Starting  ri-cell  index 

WYCEL2 

Ending  ri-cell  index 

WZCEL1 

Starting  sigma  layer  index 

WZCEL2 

Ending  sigma  layer  index 

WPRINT 

Printing  interval  (number  of  time  steps) 

WPRSTR 

Time  step  when  printing  starts 

WPREND 

Time  step  when  printing  ends 

WPRVAR 

Character  string  indicating  variables  printed 

B1 
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Note  :  The  following  characters  are  used  in  WPRVAR  for  designating 
different  variables. 


E  :  Surface  elevation  (cm) 

X  :  X-direction  unit  flow  rate  (cm^sec) 

Y  :  Y-direction  unit  flow  rate  (cm^/sec) 

U  :  X-direction  velocity  (cm/sec) 

V  :  Y-direction  velocity  (cm/sec) 

W  ;  Z-direction  velocity  (cm/sec) 

S  :  Salinity  (ppt) 

T  :  Temperature  (deg  C) 

A  :  Average  velocity  magnitude  (cm/sec)  and  direction 
(measured  clockwise  from  the  true  North,  deg) 

DUMMY 

SNPCRD  (918, A8) 

SNPCRD  :  Number  of  snapshot  control  lines  to  follow 
DUMMY 

SXCEL1 ,  SXCEI_2,  SYCEL1,  SYCEL2,  SZCEL1,  SZCEL2,  SNPINT, 
SNPSTR,  SNPEND,  SNPVAR  (9I8,A8) 

If  SNPCRD  >  0,  SNPCRD  lines  have  to  be  furnished  below. 

SXCEL1  ;  Starting  ^-cell  index 

SXCEL2  ;  Ending  §-cell  index 

SYCEL1  ;  Starting  T)-cell  index 

SYCEL2  ;  Ending  T]-cell  index 

SZCEL1  ;  Starting  sigma  layer  index 

SZCEI_2  ;  Ending  sigma  layer  index 

SNPINT  ;  Snapshot  interval  (number  of  time  steps) 

SNPSTR  ;  Time  step  when  snapshots  start 
SNPEND  :  Time  step  when  snapshots  end 
SNPVAR  ;  Character  string  indicating  snapshot 
variables  (same  notation  is  used  as  in 
WPRVAR) 

DUMMY 
NRANG  (9I8,A8) 

NRANG  ;  Number  of  ranges  for  computing  discharges 
DUMMY 

RANGDR,  RPOS1,  RPOS2,  RPOS3,  RRNAME  (7X,A1 ,3i8,A45) 

If  NRANG  >  0,  NRANG  lines  have  to  be  furnished  below. 

RANGDR  ;  Range  direction  (X  for  §  and  Y  for  tj) 

RPOS1  :  (ti)  cell  index  of  range  line 
RPOS2  ;  Starting  t]  (|)  cell  index  for  range 
RPOS3  ;  Ending  r|  (§)  cell  index  for  range 
RRNAME  :  Range  descriptor  (name) 
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DUMMY 

IGI,  IGH,  IGT,  IGS,  IGU,  IGW,  IGC,  IGQ,  IGP  (1018)  ;  Printout  flags.  A 

value  of  1  turns  printing  on  and  0  turns  it  off. 

IGI 

Print  arrays  such  as  NS,  MS,  NR,  MR,  etc. 

IGH 

Print  all  depth  arrays 

IGT  =  0 

IGS  =  0 

Print  restart  arrays 

IGU  =  0 

IGW=  0 

IGC 

Print  grid  coordinates  and  depths 

IGQ  =  0 

IGP 

Save  grid  information  in  FILE  23  for  plotting 
snapshots 

DUMMY 

XREF,  ZREF,  UREF,  COR,  GR,  ROO,  ROR,  TO,  TR  (10F8.0) 

XREF 

Reference  horizontal  grid  distance 
(Maximum  horizontal  dimension  divided  by 
number  of  cells  in  that  direction,  cm) 

ZREF 

Reference  depth  (average  depth  in  cm) 

UREF 

Reference  horizontal  velocity 
(average  velocity  in  cm/sec) 

COR 

Coriolis  parameter 

GR 

Gravitational  acceleration  (cm/sec^ 

ROO 

Minimum  density  expected  (gm/cc) 

ROR 

Reference  density  (maximum  expected)  (gm/cc) 

TO 

Minimum  temperature  (Celsius) 

TR 

Reference  temperature  (maximum  expected) 

(Celsius) 

DUMMY 

THETA  (10F8.0) 

THETA 

Time  level  weighting  factor  in  computations 
(A  value  of  1.0  was  used  in  the  Bight  model) 

DUMMY 

ITEMP,  ISALT,  ICC, 

IFI,  IFA,  IFB,  IFC,  IFD  (1018) 

ITEMP  =  0 

No  computation  of  temperature 

=-1 

Compute  temperature  (use  daily  equilibrium 
temperature  as  river  boundary  temperature) 

=-2 

Compute  temperature  (use  time-varying 
temperature  as  river  boundary  temperature) 

ISALT  =  0 

No  computation  of  salinity 

=-2 

Compute  salinity,  setting  salinity  and 
temperature  at  tidal  boundaries 

ICC  =  0 

Not  used 

IFI  =  1 

Compute  nonlinear  (inertia)  terms 

=  0 

No  computation  of  nonlinear  terms 

IFA  =  0 

Not  used 

IFB  =  0 

Not  used 
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IFC  =  0 

;  Not  used 

IFD  =  1 

;  Compute  horizontal  diffusion  terms 

=  0 

:  No  computation  of  horizontal  diffusion  terms 

DUMMY 

IVIS,  IWQ,  IPT  (1018) 

IVIS  =  1 

;  Save  information  for  subsequent  visualization 

=  0 

;  Information  for  visualization  is  not  saved 

IWQ  =  1 

:  Save  information  for  input  to  Water  Quality  Model 

=  0 

;  Information  for  input  to  Water  Quality  Model  is 
not  saved 

IPT  =  1 

;  Save  information  for  input  to  the  Particle 

Tracking  Model 

=  0 

I  Information  for  input  to  the  Particle  Tracking 

Model  is  not  saved 

DUMMY 

VISSTRF,  VISENDF,  VISINTF  (1018)  (to  be  supplied  only  if  IVIS=1) 

VISSTRF 

;  Starting  time  step  (integer)  for  saving 

information  for  flow  (surface  elevation  and 
horizontal  velocities) 

VISENDF 

;  Ending  time  step  (integer)  for  saving 
information  for  flow  (surface  elevation  and  horizontal 
velocities) 

VISINTF 

:  Time  step  interval  (integer)  for  saving  flow 
information 

DUMMY 

VISSTRS,  VISENDS,  VISINTS  (1018)  (to  be  supplied  only  if  IVIS=1) 

VISSTRS 

:  Starting  time  step  (integer)  for  saving 
salinity  information 

VISENDS 

I  Ending  time  step  (integer)  for  saving 
salinity  information 

VISINTS 

;  Time  step  interval  (integer)  for  saving 

salinity  information 

DUMMY 

VISSTRT,  VISENDT,  VISINTT  (1018)  (to  be  supplied  only  if  IVIS=1) 

VISSTRT 

:  Starting  time  step  (integer)  for  saving 
temperature  information 

VISENDT 

;  Ending  time  step  (integer)  for  saving 
temperature  information 

VISINTT 

;  Time  step  interval  (integer)  for  saving 

temperature  information 

Note:  The  ending  time  step  should  be  greater  than  the  starting  time  step. 

If  you  want  to  save  only  some  visualization  information  (e.g.,  flow)  but  not 
others  (e.g.,  salinity),  set  VISSTRF  and  VISENDF  to  realistic  values 
corresponding  to  the  simulation  but  set  VISTRS  and  VISENDS  to  values  > 

IT2. 
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DUMMY 

BVR,  S1,  S2,  PR,  PRV,  TWE,  TWH,  FKB,  TOO  (10F8.0) 


BVR  =  1.0 

:  Reference  turbulent  thermal  eddy  diffusivity 
(not  used) 

SI  =  10.0 

:  Empirical  constant  used  in  computation  of 
simple  variable  vertical  eddy  viscosity 

S2  =  3.33 

:  Empirical  constant  used  in  computation  of 
simple  variable  vertical  eddy  diffusivity 

PR  =  1. 

;  Turbulent  Prandtl  number 

PRV  =  1. 

:  Vertical  turbulent  Prandtl  number 

TWE 

;  Temperature  in  the  epilimnion  (for  computing 
initial  conditions) 

TWH 

;  Temperature  in  the  hypolimnion  (for  computing 
initial  conditions) 

FKB 

:  Vertical  grid  index  of  the  initial  thermocline  location 
(for  computing  initial  conditions) 

TOO  =  0.0 

;  Initial  surface  heat  flux  (cal/cm^sec) 

Note:  The  initial  conditions  computed  using  TWE,  TWH,  and  FKB  are 

overridden  by  FILE  17. 

DUMMY 

IVER,  ICON,  lUBO,  IBL,  IBR,  JBM,  JBP  (1018) 

IVER  =  2 
ICON  =  3 
IUBO=  0 
IBL  =  1 
IBR  =  ICELLS 
JBM  =  1 
JBP  =JCELLS 

DUMMY 

CREF,  CMAX,  CO  (10F8.0) 

CREF  =  1. 

CMAX=  100. 

CO  =  0. 

DUMMY 

ICC1,  ICC2,  JCC1,  JCC2,  ID1,  ID2,  JD1,  JD2  (1018) 

ICC1  =  0 
ICC2  =  0 
JCC1 =  0 
JCC2=  0 
ID1  =  0 
ID2  =  0 
JD1  =0 
JD2  =0 


Appendix  B  List  of  Input  Data  in  RIe  4 


DUMMY 

lEXP,  lAV,  AVR,  AVI,  AV2,  AVM,  AVM1,  AHR  (2I8,8F8.0) 

lEXP 

;  Vertical  eddy  coefficient  flag 

lEXP  =  0 

:  Constant  eddy  coefficient.  Set  ISPAC(9)=0. 

=-1 

:  Munk-Anderson  type  first  order  turbulence 
model.  Richardson  number  dependent  eddy 
coefficient  with  length  scale  linearly 
increasing  from  bottom  and  surface 

=-2 

;  Munk-Anderson  type  first  order  turbulence 
model.  Richardson  number  dependent  eddy 
coefficient  with  length  scale  linearly 
increasing  from  bottom  to  surface 

=-3 

;  Second-order  turbulence  model 

Note:  For  lEXP  <  0,  set  ISPAC(9)=1. 

lAV 

:  Reference  vertical  eddy  viscosity  flag. 

lAV  =  0 

;  Input  parameter  AVR  is  used  as  reference 
eddy  viscosity 

=  1 

:  Reference  vertical  eddy  viscosity  is  computed 
from  AV1+TXY*AV2,  where  TXY  is  the  total  wind 
stress  and  AV1  and  AV2  are  input  parameters 
(applies  only  for  IWIND=0  and  IWIND=1) 

AVR 

;  Reference  vertical  eddy  viscosity  (cm^/sec) 

AVI 

;  Background  vertical  eddy  viscosity  when  wind 
is  zero  (cm^/sec). 

AV2 

:  If  IAV=1 ,  unstratified  vertical  eddy 
viscosity  is  computed  from  AV1+TXY*AV2. 

AVM 

:  Minimum  allowable  vertical  eddy  viscosity 
(cm7sec) 

AVM1 

;  Minimum  allowable  vertical  eddy  diffusivity 
(cm®/sec) 

AHR 

;  Reference  horizontal  eddy  viscosity  or 
diffusivity  (cm^/sec) 

DUMMY 

FM1,  FM2,  ZTOP,  SLMIN,  QQMIN  (10F8.0) 

FM1 

;  Parameter  in  Richardson  number  dependent  eddy 
viscosity 

FM2 

;  Parameter  in  Richardson  number  dependent  eddy 
diffusivity 

ZTOP 

:  Distance  between  the  top  of  the  computational 
domain  and  the  free  surface  (cm).  Used  in 
computing  turbulence  length  scale. 

SLMIN 

:  Minimum  value  of  turbulence  macroscale  (cm) 

QQMIN 

:  Minimum  value  of  turbulent  kinetic  energy 
(gm/cm/sec^ 
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DUMMY 

ICUT,  KSMALL,  QCUT,  GAMAX,  GBMAX,  F2S  (2I8,8F8.0) 

ICUT  =  0 

Eddy  coefficients  constant  below  halocline 

=  1 

Eddy  coefficients  computed  below  halocline 

KSMALL 

Number  of  times  eddy  viscosity/diffusivity  in 
turbulence  model  are  smoothed  (e.g.,  5) 

QCUT 

Coefficient  in  second-order  turbulence  model 
(0.15  -  0.25) 

GAMAX 

Maximum  value  of  eddy  viscosity  (cm^/sec) 

GBMAX 

Maximum  value  of  eddy  diffusivity  (cm^/sec) 

FZS 

Turbulence  scale  is  not  allowed  to  exceed  the 
product  of  FZS  and  the  depth 

DUMMY 

IWIND,  TAUX,  TAUY  (I8,5F8.0) 

IWIND  =  0 

Steady  and  uniform  wind  stress 

=  1 

Steady  and  uniform  wind  speed 

=  2 

Steady  and  space  variable  wind  stress 

=  3 

Steady  and  space  variable  wind  speed 

=  4 

Time  variable  and  uniform  wind  stress 

=  5 

Time  variable  and  uniform  wind  speed 

=  6 

Time  and  space  variable  wind  stress 

=  7 

Time  and  space  variable  wind  speed 

TAUX 

Uniform  wind  stress  in  x-direction  if  IWIND  =0 

Uniform  wind  speed  in  x-direction  if  IWIND=1 

TAUY 

Uniform  wind  stress  in  y-direction  if  IWIND  =0 

Uniform  wind  speed  in  y-direction  if  1WIND=1 

DUMMY 

ISPAC(I),  1=1,10  (1018) 

ISPAC(I)  -  ISPAC(3)  =  0  (Not  used) 

ISPAC(4) 

=  1  ;  Flag  for  computing  open  boundary 

velocities 

ISPAC(5)  -  ISPAC(8)  =  0  (Not  used) 

ISPAC(9) 

=  0  :  For  lEXP  =  0 

=  1  ;  For  lEXP  <  1 

ISPAC(IO) 

=  0  (Not  used) 

DUMMY 

JSPAC(I),  1=1,10  (1018) 

JSPAC(I)  =  0  (Not  used) 

JSPAC(2) 

Flag  for  3-D  mode,  quadratic  friction 

=  0 

Constant  bottom  friction  factor  =  CTB 

=  1 

Bottom  friction  based  on  logarithmic  law 

JSPAC(3) 

Flag  for  Coriolis  terms 

=  0 

Coriolis  effects  accounted  for 

=-1 

Coriolis  effects  neglected 

JSPAC(4)  -  JSPAC(IO)  =  0  (Not  used) 
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DUMMY 

RSPAC(I),  1=1,10  {10F8.0) 

RSPAC(I)  - 

RSPAC(6)  =  0.  (Not  used) 

RSPAC(7) 

:  Depth  (cm)  below  which  the  bottom  friction 
coefficient  follows  a  ramp  function 

RSPAC(8) 

;  Dummy  parameter  (set  to  0.008  ) 

RSPAC(9)  - 

RSPAC(IO)  =  0.  (Not  used) 

DUMMY 

IBTM,  HADD,  HMIN,  HI,  H2,  SSSO,  HMAX  (I8,5F8.0) 

IBTM 

;  Bottom  bathymetry  flag 

IBTM  =  0 

;  Bottom  depth  varies  linearly  from  west  to 
east  of  the  basin 

=  1 

;  Bottom  depth  varies  linearly  from  south  to 
north  of  the  basin 

=  2 

:  Bottom  depth  array  for  cell  center  depths 
read  from  input  file  (FILE  4) 

=  3 

;  Bottom  depth  arrays  HS,  HU,  HV  read  from 

FILE  12 

=  4 

:  Bottom  depths  and  coordinates  of  cell  comers 
read  from  FILE  1 5  (set  ITRAN=2) 

HADD 

;  A  constant  depth  added  to  the  depth  array 
(cm) 

HMIN 

:  Minimum  water  depth  (cm) 

HI 

;  Bottom  depth  (cm)  along  the  west  or  south 
boundary  of  the  basin  for  IBTM  =  0  or  1 

H2 

;  Bottom  depth  (cm)  along  the  east  or  north 
boundary  of  the  basin  for  IBTM  =  0  or  1 

SSSO 

;  Initial  water  surface  elevation  (cm) 

HMAX 

:  Maximum  water  depth  (cm)  allowed 

DUMMY 

ISMALL,  ISF,  ITB,  ZREFBN,  CTB,  BZ1,  ZREFTN,  TZ1  (3I8,7F8.0) 

ISMALL  =  0 

:  Small  amplitude  assumption  is  invoked. 

Surface  elevation  is  not  added  to  the  still 
water  depth  to  obtain  the  total  depth 

=  1 

:  Small  amplitude  assumption  is  not  invoked. 

Surface  elevation  is  added  to  the  still 
water  depth  to  obtain  the  total  depth 

ISF  =  0 

:  Free  surface  current  flag  (set  to  0) 

ITB 

:  Bottom  friction  flag 

=  1 

;  Linear  bottom  friction  for  internal  mode 

>  1 

;  Quadratic  bottom  friction  for  internal  mode 

ZREFBN 

;  Reference  height  above  bottom  (cm) 

CTB 

:  Constant  bottom  drag  coefficient  (typical 
value  0.003) 

BZ1 

;  Bottom  roughness  height  (cm) 

ZREFTN 

;  Reference  height  at  the  top  (cm) 
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TZ1 

:  Constant  surface  roughness  height  (cm) 

DUMMY 

XMAP,  ALXREF,  ALYREF  (10F8.0) 

XMAP 

:  Mapping  factor  that  scales  the  (x,y) 
coordinates  created  by  the  grid  generation 
code  to  the  real  world 

ALXREF 

;  X-reference  length  in  the  computational  plane 

ALYREF 

:  Y-reference  length  in  the  computational  plane 

Note  :  ALXREF  and  ALYREF  are  used  if  ITRAN  =  0 

DUMMY 

ITRAN  (1018) 

ITRAN  =  0 

;  Cartesian  grid 

=  1 

:  Curvilinear  grid  created  by  WESCOR.  Cell 
corner  coordinates  read  from  FILE  15 

=  2 

:  Curvilinear  grid  created  by  WESCORA  or  EAGLE. 

Cell  corner  coordinates  and  depths  read  from 

FILE  15 

DUMMY 

ITBRK(I),  1=1,10  (10I8) 

ITBRK(I),  1= 

=1,10;  Time  steps  at  which  information  is  written  to 

hot-start  files  (increasing  order) 

DUMMY 

NSTA,  NFREQ 

NSTART  (1018) 

NSTA 

:  Number  of  stations  where  information  is  saved 

for  time  series  plots  of  currents 

NFREQ 

;  Time  step  interval  for  saving  currents 

NSTART 

:  Beginning  time  step  for  saving  currents 

DUMMY 

IST(K),  JST(K), 

STATID(K)  (214, A48) 

If  NSTA  >  0,  NSTA  lines  have  to  be  furnished  below. 

IST(K),  JST(K)  :  Cell  indices  (l,J)  of  a  station  where 

currents  are  saved 

STATID(K) 

:  Station  descriptor 

DUMMY 

NSTAS,  NFREQS,  NSTRTS  (1018) 

NSTAS 

:  Number  of  stations  where  water  surface 
elevations  are  saved  for  time  series  plots 

NFREQS 

:  Time  step  interval  for  saving  water  surface 

elevations 

Appendix  B  List  of  Input  Data  in  File  4 


B9 


NSTRTS 

;  Beginning  time  step  for  saving  water  surface 
elevations 

DUMMY 

ISTS(K),  JSTS(K),  STATS  (K)  {2I4.A48) 

If  NSTAS  >  0,  NSTAS  lines  have  to  be  furnished  below. 

ISTS(K),JSTS(K) 

;  Cell  indices  (l,J)  of  a  station  where  water 
surface  elevations  are  saved 

STATS(K) 

;  Station  descriptor 

DUMMY 

MSTA,  MFREQ,  MSTART  (1018) 

MSTA 

;  Number  of  stations  where  salinity  and 
temperature  information  is  saved  for  time 
series  plots 

MFREQ 

;  Time  step  interval  for  saving  information 

MSTART 

;  Beginning  time  step  for  saving  information 

DUMMY 

ISTSA(K),  JSTSA(K),  STATSA(K)  (214, A48) 

If  MSTA  >  0,  MSTA  lines  have  to  be  furnished  below. 

ISTSA(K),JSTSA(K) 

;  Cell  indices  (l,J)  of  a  station  where 

salinity  and  temperature  are  saved 

STATSA(K) 

:  Station  descriptor 

DUMMY 

NRIVER 

;  Number  of  river  boundaries  (218, F8. 0,418) 

NRIVER  =  0 

:  No  river  boundaries 

<  0 

;  River  inflows  are  steady 

>  0 

;  Time  variable  inflows 

If  NRIVER  =  0,  use  the  following  lines 

DUMMY 

DUMMY 

If  NRIVER  >  0,  use  the  following  lines 

DUMMY 

IJRDIR(K),  IJRROW(K),  IJRSTR(K).  IJREND(K)*  (1018) 

IJRDIR(K)  =  1 

;  River  boundary  is  on  left  (west) 

=  2 

;  River  boundary  is  on  bottom  (south) 

=  3 

;  River  boundary  is  on  right  (east) 

=  4 

;  River  boundary  is  on  top  (north) 
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IJRROW(K) 


;  Index  of  the  row  (J)  or  column  (I) 
of  the  river  boundary 
IJRSTR(K)  :  Starting  I  or  J  index  of  the  river  boundary 

IJREND(I^  :  Ending  I  or  J  index  of  the  river  boundary 

*  NRIVER  lines  have  to  be  furnished 


If  NRIVER  <  0,  use  the  following  lines 


DUMMY 

IJRDIR(K),  IJRROW(K),  IJRSTR{K),  IJREND(K)'  (1018) 


IJRDIR(K)  =  1 
=  2 
=  3 
=  4 

IJRROW(K) 

IJRSTR(K) 

IJREND(K) 


River  boundary  is  on  left  (west) 

River  boundary  is  on  bottom  (south) 

River  boundary  is  on  right  (east) 

River  boundary  is  on  top  (north) 

Index  of  the  row  (J)  or  column  (I) 
of  the  river  boundary 

;  Starting  I  or  J  index  of  the  river  boundary 
;  Ending  I  or  J  index  of  the  river  boundary 


*|NRIVER|  lines  have  to  be  furnished  followed  by  the  lines 
shown  below. 

DUMMY 

ICELL,  JCELL,  QRIVER(K,1J)*  (2I8,F8.0,4I8) 

ICELL,  JCELL  ;  Coordinates  of  a  cell  (l,J)  where  DRIVER  is 

prescribed 

QRIVER(K,IJ)  :  Steady  river  inflow 

‘Repeat  for  all  the  river  cells,  in  order. 


DUMMY 
NEAR  (1018) 

NEAR  ;  Number  of  interior  thin-wall  barriers 


If  NEAR  =  0,  use  the  following  line 
DUMMY 

If  NEAR  >  0,  use  the  following  lines 


DUMMY 

IJEDIR(l<g,  IJEROW(K),  IJESTR(K),  IJEEND(K)'  (1018) 


IJEDIR(K)  =  1 
=  2 

IJEROW(K) 

IJESTR(K) 

IJEEND(K) 


Earner  is  in  ^-direction 

Earner  is  in  ii-direction 

Index  of  row  (J)  or  column  (I)  of  barrier 

Starting  I  or  J  index  of  barrier 

Ending  I  or  J  index  of  barrier 


‘NEAR  lines  have  to  be  furnished. 
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DUMMY 

TIDFNO,  TIDBND  (1018) 

TIDFNO  ;  Number  of  tidal  elevation  tables  entered  as 

input 

TIDBND  ;  Number  of  tidal  elevation  boundaries 

DUMMY 

If  TIDFNO  >  0,  read  the  following  line(s) 

TIDSTR(I),  1=1  .TIDFNO  (10I8) 

TIDSTR(I)  :  The  entry  number  in  each  tidal  elevation 

table  corresponding  to  the  starting  time  of 
the  simulation 


DUMMY 

If  TIDBND  >  0,  TIDBND  lines  of  the  following  format  have  to  be  read. 
IJDIR(I),  IJROW(I),  IJSTR(I).  IJEND(I),  TIDTYP(I),  TIDFNI(I), 

TIDFN2(I)  (418, A8, 518) 

IJDIR(I)  =  1  :  Tidal  boundary  is  on  left  (west) 

=  2  ;  Tidal  boundary  is  on  bottom  (south) 

=  3  :  Tidal  boundary  is  on  right  (east) 

=  4  ;  Tidal  boundary  is  on  top  (north) 

IJROW(I)  ;  Index  of  the  row  (J)  or  column  (I)  of  the 

tidal  boundary 

IJSTR(I)  :  Starting  I  or  J  index  of  the  tidal  boundary 

IJEND(I)  ;  Ending  I  or  J  index  of  the  tidal  boundary 

TIDTYP(I)  =  "CONSTANT'  ;  Constant  tidal  elevation  between 

IJSTR(I)  and  IJEND(I) 

=  "INTERP"  ;  Linear  interpolation  of  tidal 

elevation  between  IJSTR(I)  and 
IJEND(I) 


TIDFNI(I) 

TIDFN2(I) 


;  The  number  of  the  tidal  elevation  table  for 
CONSTANT  or  INTERP  type  of  boundaries 
;  The  number  of  the  second  tidal  elevation 
table  used  for  interpolation  on  INTERP  type 
boundaries 


Optional  input: 


DUMMY 
l.J  (4x, 12, lx, 12) 


;  Indices  of  a  cell  where  HS  is  reset  to  0. 


DUMMY 
1,J  (4x,l2,1x,l2) 

DUMMY 
1,J  (4x,12,1x,l2) 


;  Indices  of  a  cell  where  HU  is  reset  to  0. 


;  Indices  of  a  cell  where  HV  is  reset  to  0. 


DUMMY 

1,J,  RDEPTH  (free  format)  ;  Indices  and  depth  (ft)  of  a  cell  where 

HS  is  reset  to  non-zero  value  RDEPTH. 
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List  of  Input  Data  Files^ 


FILE  13 

River  inflows  are  read  from  FILE  13.  These  data  are  read  first  as  a 
time  line  {DAY  and  HOUR)  formatted  by  218.  Next,  the  (l,J)  location  and 
discharge  in  cubic  feet  per  second  for  each  cell  of  each  river  boundary 
are  read  and  formatted  by  (218,  F8.0). 

FILE  14 

Wind  data  are  read  from  FILE  14.  These  data  are  in  the  form  of  time 
(DAY  and  HOUR)  and  the  x  and  y  components  of  the  wind  velocity  in 
meters  per  second  of  each  wind  field  used.  These  data  are  formatted  by 
(2I5,6F10.0). 

FILE  15 

The  (x,y)  coordinates  and  depths  of  the  New  York  Bight  grid  cell 
corners  are  read  from  FILE  15.  This  file  was  created  from  a  run  of  the 
grid  generation  code  EAGLE  and  a  depth  interpolation  program.  The  first 
line  contains  the  file  name  formatted  as  A80.  The  number  of  corner 
points  in  f  and  /?  are  read  next  unformatted.  The  coordinates  and 
depths  are  read  next  unformatted,  one  line  per  corner. 

FILE  16 

Tabular  tide  data  are  read  from  FILE  16.  The  first  line  is  the  title 
formatted  as  A80.  The  tide  data  are  in  the  form  of  time  (MONTH,  DAY, 
YEAR,  HOUR,  MINUTES)  and  the  water  surface  elevations  in  centimeters 
relative  to  selected  datum  for  TIDFNO  points.  These  data  are  formatted 
by  (I2,1X,2I3,1X,2I2,(T17,8F8.2)). 


'  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  (metric)  units  is 
presented  on  page  vi. 
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FILE  17 

The  initial  temperature  field  in  degrees  Celsius  is  read  from  FILE  17 
by  format  (1 OE1 2.5).  This  file  is  created  from  a  few  observed  values.  The 
resulting  field  is  then  smoothed  in  the  f  and  rj  directions  several  times  before 
it  is  written  to  FILE  17. 

FILE  18 

The  Manning's  roughness  n  is  hard-wired  in  the  model  for  the  New 
York  Bight  application.  By  changing  the  source  code,  a  field  of  Manning's  n 
values  may  be  input  by  format  (20F4.0).  The  input  values  are  multiplied  by 
0.001  in  the  source  code  to  yield  the  actual  values.  They  are  input  by  rows. 

FILE  19 

Daily  average  equilibrium  temperatures  in  degrees  Celsius  and  surface 
heat  exchange  coefficients  in  units  of  cm/sec  are  read  from  FILE  19.  These 
data  are  in  the  form  of  time  (DAY  and  HOUR),  equilibrium  temperature,  and 
heat  exchange  coefficient.  They  are  formatted  by  {2I5,F1 0.0,E1 2.5). 

FILE  74 

The  initial  salinity  field  in  parts  per  thousand  is  read  from  FILE  74  by 
format  (10E12.5).  This  file  is  created  in  the  same  fashion  as  FILE  17. 

FILE  76 

Time-varying  salinity  in  parts  per  thousand  and  temperature  in 
degrees  Celsius  at  tidal  boundaries  are  read  from  FILE  76  if  salinity  and 
temperature  are  to  be  computed.  These  data  are  in  the  form  of  time  (DAY 
and  HOUR)  formatted  by  (215).  Next,  the  (l,J)  location  of  each  tidal 
boundary  cell  and  the  vertical  distribution  of  salinity,  starting  from  the  top 
layer  to  the  bottom  layer  are  read.  These  data  are  followed  by  temperature 
data  using  the  same  format  as  for  the  salinity.  The  format  is  (215,1 1  F5.0) 

FILE  78 

Time-varying  temperature  data  at  river  flow  boundaries  are  read  from 
FILE  78  if  temperatures  are  to  be  computed  and  equilibrium  temperatures  are 
not  used  as  river  boundary  temperatures.  These  data  are  in  the  form  of  a 
time  (DAY  and  HOUR)  formatted  by  (215).  Next,  the  (l,J)  location  of  river 
flow  boundary  cells  and  corresponding  temperatures  starting  from  top  layer 
to  bottom  layer  are  read.  These  data  are  formatted  by  (215,1 1  F6.0). 

FILE  85 

This  file  contains  corner  point  information  for  the  water  quality  boxes 
(free  formatted). 
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FILE  90 


This  one  line  file  contains  the  time  averaging  interval  (NAVG)  and  the 
starting  time  step  (ITPTS)  for  saving  information  for  input  to  the  Particle 
Tracking  Model  (free  formatted). 

FILE  95 

This  file  contains  the  number  of  surface  boxes  (NSB),  the  time 
averaging  interval  (NAVG),  the  time  step  to  start  saving  information 
(ITWQS),  WQMINT,  WQMSTR,  WQMEND,  and  WQMSNP  (free  format).  It 
also  has  information  on  each  of  the  surface  boxes  and  each  flow  face. 
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Appendix  D 

Input  Data  in  File  4  for 

April  1976  Application 


TITLE(A80) 

NY  BIGHT  TEST 

:  76  X  4S  GRID  RERUN  18 

IFD  = 

1?  3-D  »  Ue  SOdays  B-13-93 

ITl 

IT2 

DT 

ISTtffT 

ITEST 

nSALT 

1 

17280 

ISO.O 

0 

0 

2880 

HPRCHD 

NXCEU 

inxa2 

NYCELl 

KYCEL2 

WCELl 

WCEL2 

NPRIHT 

NPRSTR 

KPREND  IFRVAR 

1 

76 

1 

45 

1 

10 

17280 

72800 

100000  EST 

SNPcra) 

i 

l 

sum 

SXCEL2 

SYCEU 

SYCEL2 

SZCELl 

SICEL2 

SW>IKT 

SIfSTR 

SNPVAR 

1 

76 

1 

45 

1 

10 

17280 

72800 

100000  EST 

M»NG 

RANGOR 

RPOSl 

fiP0S2 

RP0S3 

RRIWC 

X 

40 

40 

40 

EAST  RIVER(HEST) 

1 

44 

40 

40 

EAST  RIVERIEASn 

IGI 

lat 

IGT 

I6S 

IGU 

IGM 

IGC 

IGQ 

IGP 

0 

0 

0 

0 

0 

0 

0 

0 

0 

XREF 

ZKEF 

UREF 

COR 

GR 

ROO 

ROR 

TO 

TR 

S21S00. 

9144. 

10.0 

.00009 

981.0 

1.0 

1.021 

1. 

35. 

■noA 

1.0 

ITEJP 

ISALT 

ICC 

IFI 

IFA 

IFB 

IFC 

IFD 

AHO 

-2 

-2 

0 

1 

0 

0 

0 

1  100000. 

IVIS 

!HQ 

IPT 

0 

1 

0 

VISSTRF  VISENDF 

VISINTF 

VISSTBS  VISENDS  VISIKTS 

VISS7RT  VISENDT  VISINH 

GVR 

SI 

S2 

PR 

PRV 

WE 

TW 

FO 

TOO 

1.0 

10. 

3.33 

1. 

1.0 

10. 

10. 

S. 

0.0 

IVER 

ICON 

lueo 

IfiL 

I6R 

JEH 

3BP 

2 

3 

0 

1 

76 

1 

45 

CREF 

OKI 

CO 

1. 

100. 

0. 

IICl 

IC2 

JJCl 

JC2 

IDl 

ID2 

JDl 

JD2 

0 

0 

0 

0 

0 

0 

0 

0 

Appendix  D  Input  Data  in  File  4  for  April  1 976  Application 


D1 


•J 

V 

lU. 

FHl 

FH2 

2TEP 

-.5 

-1.5 

0. 

lOIT 

XSHMl 

OOIT 

0 

0 

0.15 

[HIND 

TAUX 

TAUY 

5 

0.00 

0.00 

0. 

0. 

10.0 

SLNIN 

QGHIN 

1. 

0.01 

F2S 

5000. 

5000. 

0.2 

0.005  100000. 


ISPAC(I),M.10 
0  0 

JSPfiCCD.M.lO 
1  1 

RSPAC(I).1=1.10 


.020 

0. 

.00001 

.00001 

0. 

1. 

1000. 

IBTH 

J 

■sa 

HHIN 

HI 

H2 

(flAX 

4 

u. 

304.8 

0.0 

0.0 

0.0  999999. 

isim. 

ISF 

ITB 

ZREFBN 

CTB 

BZl 

ZREFBN 

1 

0 

5 

5. 

0.0025 

0.005 

5. 

XMAP 

ALXREF 

ALYREF 

30.48 

521500 

747400 

ITRAN 

IBD(l) 

(2) 

(3) 

(4) 

2 

4 

4 

4 

4 

ITBRXd) 

(2) 

(3) 

(4) 

(5) 

(6) 

70272  88128  105408  0  0 

NSTA  NFREQ  KSTART  (CURRENT  STATIONS) 
8  48  0 

1ST  JST  STATID(K)  (2I4.A48)  (  ONE  CAR 

37  15  lESA  LT  1 

7  15  HESALT2 
11  6  HESA  LT  3 

54  17  fESALT4 
52  9  BESALT5 
23  24  lESA  LT  6 
40  42  HUDSON  R.  (40.42) 

38  41  HUDSON  R.  (38.41) 

NSTAS  IfREQS  NSTRTS  (TIDE  STATIONS) 


(  ONE  CARD  FDR  EACH  STATION  ) 


13  48  0 

1ST  jgr  STAnD(K)  (214, A4a) 

4  25  ATUWnC  cm  (4,25) 

27  37  SAfOY  HOOK  (27,37) 

39  41  THE  BATTERY  (39,41) 

62  37  JEN  LONDON  (62.37) 

68  43  NENFORT  (68,43) 

40  42  HUDSON  R.  (40,42) 

38  41  HUDSON  R,  (38,41) 

,44  41  L.I.  SOUND  (44,41) 

33  7  OCEAN  BDRY  (33.7) 

75  33  (75,33) 

74  33  (74,33) 

75  35  (75,35) 

75  31  (75,31) 

HSTA  KTEQ  KST^TT  (SAL 

7  144  2880 

1ST  JST  STATID(K)  (2I4.A48) 
37  15  HESA  LT  1  (37,15) 

7  15  lESA  LT  2  (7,15) 

11  6  HESA  LT  3  (11,6) 


(  0(€  CARD  FOR  EACH  STATION  ) 


(SALINITY  STATIONS) 

48)  (  Of£  CARO  FOR  EACH  STATION  > 


Appendix  D  Input  Data  in  File  4  for  April  1976  Applications 


54  17  HESALT 
52  9  »€SA  LT 

23  24  reSALT 
2&  IS  KSA  LT 
»RIVER 
1 

IJRDIR  IJRROH 
3  40 

1  3 
NBAR 

0 

IJROIR  IJRRCM 
TIERO  TIDBND 
11  8 
TIDSTR  2 
1  1 
1 

IJDIR  IJraU 
1  1 
1  1 
1  1 

2  1 

2  1 
2  1 
2  1 
3  76 

RESET  KS(I,3)  TO 
24.45 
61.26 
62.26 
63.27 

58.30 

59.30 

66.30 

59.31 

60.32 

52.33 

26.36 

27.36 

71.36 

28.37 

34.37 

35.37 

64.37 

71.37 

72.37 

62.38 

72.38 

35.39 

44.39 

45.39 

46.39 

47.39 

72.39 

72.40 

40.41 

41.41 


4  (54,17) 

5  (52.9) 

6  (23,24) 

7  (26,15) 


I^ISTR  IJREKD  (  {*£  (^  FDR  EACH  RIVER  ) 
42  42 

{RIVER  {  ONE  C;«D  FOR  EACH  CELL  ) 

IJRSTR  IJREND  (  ONE  CARD  FOR  EACH  BAR  ) 


I3STRT  I3END  TIDIYP  TIDFHl  TII1FN2 
1  6IKTERP  4  3 

4  SIKTERP  3  2 

8  25IKTERP  2  1 

1  lOIKTERP  4  5 

10  32I)(TERP  5  6 

32  60I)nERP  6  7 

60  76IKTERP  7  7 

1  43I)(TERP  7  11 

ZERO  AT  THE  FOUDHINS  CELLS 


ppendix  D  Input  Data  in  File  4  for  April  1 976  Application 


42.41 

43.41 

33.42 

36.42 

41.42 

42.42 

43.42 

49.43 

50.43 

51. 43 

52.43 

53.43 

54.43 

56.43 

reset  hj(i.j)  to  zero  at  tie  following  miQ 
61,26 
62,26 

59.30 

26.36 

27.36 

35.37 

36.37 

72.37 

63.30 

36.39 

44.39 

45.39 

46.39 

47.37 

41.41 

42.41 

43.41 

36.42 

42.42 

43.42 

50.43 

51.43 

52.43 

53.43 

54.43 

57.43 

RESET  HVd.J)  TO  ZERO  AT  THE  FOLLOWING  mig 

59.31 
52.33 

71.37 

64.38 

72.38 
.  44,39 

45.39 

46.39 

47.39 
62.37 

72.39 

72.40 

35.43 

36.43 

49.44 

50.44 


Appendix  E  River  Inflows  in  Fiie  1 3 


0 

0 

-4.745 

0.080 

-5.165 

5.655 

-5.550 

4.175 

0 

& 

-5.825 

3.270 

-6.425 

7.505 

-4,215 

7,415 

0 

12 

-6.375 

9.670 

-7.420 

9.255 

-0.770 

10,425 

0 

18 

-2.690 

10.980 

2.320 

4.315 

6.340 

4.720 

1 

0 

-0.210 

6.150 

3.410 

1.620 

3.645 

1.340 

1 

6 

2.990 

0.330 

4.905 

1.215 

4.765 

0.660 

1 

12 

2.860 

-1.620 

3.035 

-1.165 

3.425 

-1.765 

1 

18 

2.440 

-0.615 

3.520 

-1.640 

4.515 

-1.670 

2 

0 

1.875 

-2.640 

3.325 

-6.^ 

5.200 

-6.595 

2 

6 

5.800 

-9.530 

6.450 

-10.430 

7,610 

-12.375 

2 

12 

10.130 

-11.370 

9.B30 

-11.270 

10.640 

-12.920 

2 

18 

12.450 

-10.640 

11,270 

-9.200 

12.295 

-8.345 

3 

0 

10.485 

-9.780 

9,920 

-7.315 

12.700 

-4,390 

3 

6 

10,805 

-9.065 

7.970 

-4. 480 

8.525 

-6.940 

3 

12 

10.840 

-4.395 

4.210 

-1.115 

6.645 

2.870 

3 

18 

6.820 

-0.280 

0.280 

-0.275 

5.145 

5.945 

4 

0 

4.125 

1.540 

0.790 

-5,300 

2.940 

1.825 

4 

6 

7.315 

-5.945 

5.585 

-11.920 

4.640 

-8.580 

4 

12 

11.105 

-8.605 

11.070 

-10.585 

6.700 

-12.160 

4 

18 

14.455 

-6.010 

12.910 

-5.980 

9.510 

-7.545 

5 

0 

12.695 

1.205 

10.825 

I.5S5 

8.750 

1.455 

5 

6 

10.175 

1.835 

8.500 

2.470 

6.440 

3.250 

5 

12 

8.055 

2.210 

7.205 

2.130 

6.400 

3.445 

5 

18 

7.050 

4.465 

6.140 

2.950 

4.610 

3.895 

6 

0 

7.000 

4.500 

5.485 

3.520 

3.955 

4.180 

6 

6 

4.215 

0.020 

3.290 

-2.485 

2.755 

0.260 

6 

12 

4.350 

-0.650 

3.020 

-4.515 

0.510 

-3,050 

6 

18 

4.050 

-4.325 

3.695 

-5.050 

0.990 

-4.595 

7 

0 

4.250 

-2.655 

3.135 

-1.865 

-0.315 

-2.605 

7 

b 

0.685 

-3.330 

-0.355 

-2.640 

-0.770 

-2.245 

7 

12 

-0.210 

-0.320 

-0.795 

-1.735 

-2.340 

-0.520 

7 

18 

-0.920 

1.220 

-1.110 

-1.240 

-3.510 

0.555 

8 

0 

-3.490 

1.240 

-3.565 

-1.935 

-5.885 

-2.110 

8 

6 

-6.495 

-4.030 

-4.490 

-6.110 

-7.495 

-6.630 

8 

12 

-10.660 

-7.735 

-9.020 

-10.890 

-10.200 

-13.940 

8 

18 

-11.420 

-9.410 

-8.025 

-11.670 

-6.065 

-17.790 

7 

12 

2.725 

-18.945 

1.525 

-15.975 

5.470 

-15.945 

9 

18 

9.625 

-15.805 

9.840 

-10.315 

9.725 

-12,490 

10 

0 

13.495 

-2.725 

11.695 

-2,285 

8.610 

-7.430 

10 

6 

15.765 

0.410 

12.635 

1.810 

13.730 

2.500 

11 

0 

4.850 

-4.080 

5.505 

-9,045 

6.545 

0,125 

11 

b 

5.220 

-12.870 

4.600 

-13.785 

1.145 

-10.515 

11 

12 

5.695 

-13.495 

5.825 

-12.990 

4.335 

-14.300 

11 

18 

8.415 

-10.445 

9.600 

-10.380 

7.875 

-10.510 

12 

0 

8.495 

-6.250 

9.775 

-6,485 

8.080 

-5.585 

Appendix  G 

Cell  Corner  Coordinates  and 
Depths  in  File  1 5^ 


KVB  74X45  6RID  COORDINATES  AND  Cai  COW®  DEPTTIS  IN  FECT 

77  44 

485440.00 

84328.00 

2871.58 

712180.00 

134440.00 

1102.72 

739020.00 

175750.00 

845.24 

748300.00 

209820.00 

1044.44 

794390.00 

239720.00 

811.28 

819210.00 

242880.00 

494.21 

638940.00 

284540.00 

922.45 

854720.00 

302050.00 

1050.45 

870340.00 

317240.00 

1747.19 

883540.00 

328250,00 

2042.94 

894500.00 

338390.00 

1307.50 

907120.00 

348440.00 

1444.82 

920590.00 

359190,00 

1394,77 

934740.00 

370350.00 

2398.77 

949350.00 

381240.00 

1430.39 

944000.00 

391580.00 

2198.52 

979940.00 

.  400390,00 

1998.72 

993410.00 

409090.00 

•1401,43 

1007400.00 

417370.00 

2024.34 

1019400.00 

425990.00 

1485.39 

1030700.00 

433970.00 

1898.15 

1040500.00 

441280.00 

1442.43 

1049300.00 

447440,00 

1443.20 

1054500.00 

453770.00 

3408.01 

1042200.00 

459810.00 

2390.14 

1044700.00 

445400.00 

2530.58 

1070900.00 

471180.00 

3874.83 

1075700.00 

474150.00 

3357.38 

1080400.00 

480340.00 

3540.52 

1085400.00 

485320.00 

2494,58 

1090700.00 

490050.00 

3439,25 

1095900.00 

494570.00 

4432.44 

1100200.00 

499590.00 

3921.73 

1105400.00 

504250.00 

2999.45 

1109400.00 

509590.00 

2482.54 

1114300.00 

514710.00 

1953,10 

1119200.00 

519800.00 

2142.94 

1124300.00 

524940,00 

2447,65 

1130300.00 

530120.00 

2342.08 

1135900.00 

534300.00 

2228.55 

‘  A  table  of  factors  for  converting  non-Sl  units  of  measurement  to  SI  (metric)  units  is 
presented  on  page  vi. 


Appendix  G  Cell  Corner  Coordinates  and  Depths  in  RIe  15 


Appendix  H 
Tabular  Tide  in  File  16 


New  York  Bight 

B  -  Constituents  Tide 

4/01/76  -  9/30/76 

4 

1  76 

0  0 

0.00  0.00  0.00 

0.00  0.00  0.00 

0.00  0.00 

0.00 

0.00 

0.00 

4 

1  76 

030 

-39.56  -34.41  -33.85 
-26.85  -26.86  -31.15 

-32.72  -33,01 

-30.86 

-28.67 

-26.86 

4 

1  76 

I  0 

-49.44  -42.68  -41.95 
-35.19  -35.28  -41.09 

-40.35  -40.79 

-38.45 

-36.15 

-35.05 

4 

1  76 

I3C' 

-55.99  -48.03  -47.19 
-41,12  -41.30  -48.22 

-45.25  -45.79 

-43,40 

-41.12 

-40.85 

4 

1  76 

2  0 

-58.83  -50.14  -49.25 
-44.30  -44.55  -52.12 

-47.13  -47.74 

-45.43 

-43.30 

-43.92 

4 

I  76 

230 

-57.82  -48.94  -46.06 
-44.57  -44.37  -52.58 

-45.92  -46.57 

-44.46 

-42.60 

-44,10 

4 

1  76 

3  0 

-53.04  -44.53  -43.74 
-41.95  -42.26  -49.59 

-41.74  -42.40 

-40.60 

-39.10 

-41.42 

4 

1  76 

330 

-44.90  -37.26  -36,61 
-36.66  -37.00  -43.43 

-34.90  -K.54 

-34.15 

-33.09 

-36.10 

4 

1  76 

4  0 

-33.92  -27.65  -27.19 
-29.09  -2'9.42  -34.53 

-25.91  -26.49 

-25.59 

-25.01 

-2S.53 

4 

1  76 

430 

-20.94  -16.37  -16.13 
-19.77  -20.08  -23,50 

-15.38  -15.89 

-15.52 

-15.43 

-19.26 

4 

1  76 

5  0 

-6.80  -4.19  -4.20 
-9.36  -9.61  -11.09 

-4,05  -4.48 

-4.66 

-5.03 

-8.92 

4 

1  76 

530 

7.56  8.05  7.79 
1.43  1.24  l.BO 

7.32  6.97 

6.26 

5.47 

1.77 

4 

1  76 

6  0 

21.12  19.50  18.99 
11.84  11.73  14.33 

17.92  17.65 

16.46 

15.32 

12.06 

4 

1  76 

630 

33.0!  29.41  28.69 
21.21  21.17  25.66 

27.08  26.88 

25.28 

23.89 

21.29 

4 

1  76 

7  0 

42.43  37.09  36.21 
23.87  28.91  35.02 

34.16  34.00 

32.09 

30.56 

28.81 

4 

1  76 

730 

43.71  42.03  41.03 
34.31  34.43  41.79 

38.68  38.54 

36.43 

34.88 

34.13 

4 

1  76 

8  0 

51.45  43,88  42.83 
37.17  37.36  45.51 

40.33  40.19 

37.99 

36.54 

36.87 

4 

1  76 

830 

50,51  42.50  41.48 
37,24  37.48  45.90 

33.99  38.83 

36.67 

35.43 

36.85 

4 

1  76 

9  0 

45.87  37.99  37.06 
34.54  34.82  43.01 

34.76  34.55 

32.56 

31.61 

34.08 

4 

1  76 

930 

37,92  30.64  29,86 
29.24  29.53  36.97 

27.91  27.63 

25.91 

25.33 

28.75 

4 

1  76  10  0 

27.16  20.94  20.37 
21.70  21.98  28.19 

18.90  16.55 

17.18 

17.01 

21.21 

4 

1  76  1030 

14.33  9.52  9,20 
12,43  12,68  17.31 

8.34  7.90 

6.94 

7.21 

11.97 

Appendix  H  Tabular  Tide  in  Rle  16 


Appendix  I 

Initial  Temperature  Field 
in  File  17 


O.OOOOOE+OO  0.im)0EtO2  0.12500e*02  0.12500e*02  0.12500£«02  0.12500e<02  O.lZMEiOl  0.12300e402  0.125(X€«»2  0.1250C€«)2 

0.13SX£t02  0.12400E'K)2  0.12400£*02  0.12«)0E*O2  0.12600E*«2  0.12700e«)2  0.1270C£«)2  0.12700e«)2  0.12700E4<12  0.12700E«)2 

0.127«C«02  0.1Z700E«02  0.12?OC£H>2  0.12600£«fl2  0.1260Ce402  0.12600e-t02  0.12300e«02  0.123006<02  0.11?0Cet02  0.«C0Cie*Ol 

O.g3(»0E<«l  0.6«)0ce«01  O.raOOOEtCl  0.n200£«<l2  0.I23O0E«O2  0.12300E«02  0.11900£<<)2  0.12300e«02  0.12600E«Q  0,I270C£*02 

t.l2Ml£*92  t.l260Ce-K)2  0.12i0CC«<l2  0.12i>00E«<l2  O.I240CE<«2  0.12400e'H)2  0.12600E4<)2  0.12&(Xe<02  O.12iOC£402  O.12iOOE402 

«.I2590C««I2  4.1240CE402  0.12M0e«r2  0.12400E«E  0.12600E4<)2  0.1230C€4«2  0.12600E4C2  0.12700e402  0,127006402  0.124006402 

0.125006402  0.123006402  0.121006402  0.  fl9006402  0.116006402  0.117006402  0.117006402  0.118006402  0.121006402  0.123006402 

0.1230C6402  0.123006402  0.I200C6402  0.119006402  0.120006402  0.122006402  0.12000E402  0.000006400  0.O0OOC64OO  0.12S006402 

0.123006402  0.125066402  0.125006402  0.125006402  0.123006402  0.125006402  0.125006402  0.125006402  0.123006402  0,125006402 

0.124006402  0,126006402  0,126006402  0.126006402  0.127006402  0.127006402  0.127006402  0.127006402  0,127006402  0,127006402 
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0.11866£402  0.118486402  0.118646402  0.118986*02  0.118866*02  0.118496*02  0.117956402  0.117326*02  0.117036*02  0.117326*02 
0.117816402  0.118386402  0.119206*02  0.120576402  0.121926*02  0.122786*02  0.1230*402  0.122656*02  0.121666402  0.120396*02 
0.119046*02  0.117966402  0.117246*02  0.117206*02  0.117876*02  0,119046*02  0.119796402  0.119316*02  0.118116402  0.115506*02 
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Appendix  K 
Time-Varying  Vertical 
Distributions  of  Salinity 
and  Temperature  at  the  Ocean 
Boundary  In  File  76 


I  1  32.7  33.6  34.1  34.4  34.7  34.9  35.1  35.2  35.2  35.3 

I  1  16.9  15.8  15.2  15.0  14.6  14.3  13.8  13.4  13.0  12.7 

1  2  32.7  33.6  34.1  34.4  34.7  24.9  35.1  35.2  35.2  35.3 

1  2  16.9  15.8  15.2  15.0  14.6  14.3  13.3  13.4  13.0  12.7 

1  3  32.7  33.6  34.1  34.4  34.7  34.9  35.1  35.2  35.2  35.3 

I  3  16.9  15.8  15.2  15.0  14.8  14.3  13.6  13.4  13.0  12.7 

I  4  32.4  32.8  33.2  33.5  33.7  33.9  34.0  34.2  34.3  34.3 

1  4  16.7  15.9  15.4  14.8  13.3  12.0  11.0  10.6  10.4  10.3 

1  5  32.3  32.5  32.9  33.2  33.4  33.5  33.6  33.8  34.0  34.0 

1  5  16.6  15.9  15.4  14.7  12. S  11.2  10.0  9.6  9.5  9.4 

1  6  32.2  32.3  32.6  32.9  33.1  33.2  33.3  33.5  33.7  33.7 

I  6  16.5  16.0  15.5  14.6  12.5  10.6  9.3  9.0  8.9  8.8 

1  6  32.2  32.3  32.6  32.9  33.1  33.2  32.3  33.5  33.7  33.7 

1  6  16.5  16.0  15.5  14.6  12.5  10.6  9.3  9.0  8.9  8.8 

1  7  32.!  32.2  32.5  32.7  32.9  33.0  32.1  33.3  33.5  33.5 

1  7  16.5  16.0  15.5  14.6  12.4  10.5  9.1  8.6  8. 5  8.4 

1  a  32.0  32.1  32.3  32.6  32.7  32.8  32.9  33.1  33.3  33.4 

1  B  16.5  16.1  15.6  14.7  12.5  10.7  9.3  8.7  8.4  8.2 

I  a  32.0  32.1  32.3  32.6  32.7  32.8  32.9  33.1  33.3  33.4 

I  8  16.5  16.1  15.6  14.7  12.5  10.7  9.3  8.7  3.4  .8.2 

1  9  32.0  32.0  32.2  32.4  32.5  E.7  32. S  33.0  33.2  33.3 

1  9  16.5  16.1  15.7  14.8  12.9  11.2  9.8  9.0  8.6  8.1 

I  10  32.0  32.0  32.1  32.3  32.4  32.5  32.7  32.9  33.0  33.2 

1  10  16.4  16.2  15. H  15.0  13.4  11.9  10.6  9.6  9.0  8.3 

1  11  32.0  32.0  32.1  32.2  32.3  32.4  32.6  32.7  32.?  33.1 

1  II  16.4  16.2  16.0  15.2  14.0  12. B  11.5  10.2  9.5  8.6 

1  12  32.0  32.0  32.0  32.2  32.2  32.3  32.4  32.6  32.7  32.9 

1  12  16.4  16.3  16.1  15.6  14.7  13.7  12.5  !!.4  •<).:  9.0 

1  13  32.0  32.0  32i0  32.1  32.1  32.2  32.3  32.5  32.6  32.8 

1  13  16.5  16.4  16.3  15.9  15.3  14.4  13.2  12.2  10.9  9.6 

1  14  31.9  31.9  31.9  32.0  32.0  32.2  32.3  32.4  32.5  32.7 

1  14  16.6  16.6  16.5  16.3  15.6  14.9  13.7  12.S  11.6  10.4 

I  15  31.9  31.9  31.9  31.9  31.9  32.1  32.3  32.4  32.5  32.6 

1  IS  16.6  16.B  16.8  16.6  16.3  15.2  13.9  13.2  12.2  11.1 

1  16  31.8  31.8  31.8  31.9  31.9  32.1  32.3  32.4  32.4  32.5 

1  16  17.0  17.1  17.0  16.9  16.6  15.3  13.9  13.4  12.6  11.8 

1  17  31.8  31.8  31.8  31.8  31-8  32.1  32.3  32.4  32.4  32.5 

1  17  17.2  17.2  17.2  17.1  16.9  15.4  13.9  13.6  13.0  12.4 

I  IB  31.8  31.8  31. B  31.8  31. 8  32.1  32.3  12.3  32.3  32.4 

1  18  17.3  17.3  17.3  17.2  17.0  15.5  14.1  13.9  13.5  13.0 

1  19  31.7  31.7  31.7  31.8  31.8  32.0  32.2  32.2  32.2  32.3 

I  19  17.3  17.3  17.3  17.2  17.0  15.7  14.5  14.3  14.1  13.7 

1  20  31.8  31.8  31.8  31.8  31.8  32.0  32.1  32.1  32.1  32.2 

1  20  17.3  17.3  17.2  17.1  17.0  16.0  15.0  14.9  14.7  14.5 

1  21  31.8  31.8  31.8  31.8  31.8  31.9  32.0  32.0  32.0  37.1 

1  21  17.3  17.2  17.2  17.0  16.9  16.2  15.5  15.4  13.3  15.1 

I  22  31.8  31.8  31. e  31. B  31.8  31.9  31.9  31.9  31.9  32.0 

1  22  17.3  17.2  17.1  16.9  16.9  16.4  16.0  1S.B  13.8  15.6 

1  23  31.8  31.8  31.8  31.8  31.8  31.8  31.9  31.9  31.9  32.0 

1  23  17.2  17.1  17.1  16.9  16.8  16.3  16.2  16.1  16.1  16.0 

1  24  31.6  31.8  31.6  31.8  31.8  31.8  31.8  31.8  31.8  31.9 

1  24  17.2  17.1  17.0  16.8  16.8  16.6  16.4  16.3  16.3  16.2 

1  Z  31.6  31.1  31.1  31.6  31.8  31.8  31.8  31.8  31.8  31.9 

1  Z  17.2  17.1  17.0  16.8  16.8  16.7  16.6  16.S  16.3  16.4 

1  1  32.7  33.6  34.1  34.4  34.7  34.9  33.1  35.2  35.2  33.3 

1  1  16.V  15.8  13.2  15.0  14.8  14.3  13.B  13.4  1X0  12.7 

2  1  32.8  33.5  34.0  34.3  34.6  34.8  34.9  35.1  35.2  Z.2 

2  1  17.0  16.0  0.3  15.1  15.0  14.8  14.2  13.7  1X3  1X0 

3  1  32.8  3X3  34.0  34.3  34.8  34.1  34.9  35.1  3X2  3X2 
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Appendix  L 
Time- Varying  Vertical 
Distribution  of 
Temperature  at  the  River 
Boundary  in  File  78 
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